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ABSTRACT 


ADIN'A,  a  finite  element  program  for  automatic  dynamic  incremental  nonlinear 
analysis,  is  used  to  conduct  a  stress  analysis  of  the  Tarawa  class  (LHA-1)  superheater 
header.  This  investigation  is  sponsored  by  the  Naval  Sea  Systems  Command 
(NAVSEA)  as  a  result  of  cracking  in  the  superheater  header  tube  attachment  weld.  A 

2- dimensicnal  solid  axisymmetric  finite  element  model  was  developed  using  both  a 
coarse  and  fine  mesh.  Both  oressurc  and  thermally  induced  stresses  were  studied.  A 

3- dimensional  solid  finite  element  model  was  developed  and  a  preliminary  stress  analysis 


conducted. 
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I.  INTRODUCTION 


The  Tarawa  (LHA-1)  class  amphibious  assault  ships  have  been  experiencing  nu¬ 
merous  failures  of  their  superheater  header.  The  failures  occur  in  the  tube  attachment 
weld.  This  is  a  relatively  new  class  of  ship  and  they  are  expected  to  be  in  commission 
for  another  twenty  to  thirty  years.  The  solution  of  this  problem  is  imperative  to  prevent 
the  excessive  down  time  experienced  by  these  engineering  plants. 

A  2-dimensional  solid  axisymmetric  finite  element  model  was  developed  using 
ADINA.  A  study  of  the  effect  of  different  thermal  gradients  across  the  thickness  of  the 
header  was  conducted.  The  model  yielded  very  different  results  dependent  on  the 
boundary  conditions  imposed.  A  discussion  is  made  of  the  boundary  conditions  and  the 
conclusion  that  a  3-dimensional  model  would  more  closely  model  the  actual  physical 
boundary  conditions  of  this  structure. 

A  3-dimensionai  solid  finite  element  model  was  developed.  A  discussion  is  made  of 
the  boundary  conditions  and  a  study  is  made  of  the  internal  pressure  loading. 


II.  DESCRIPTION  OF  THE  PROBLEM 


NAVSEA  has  been  conducting  an  extensive  investigation  to  determine  the  cause  for 
the  weld  joint  cracking.  The  following  list  was  provided  by  NAVSEA  and  is  an  analysis 
of  weld  joint  cracking  on  the  L'SS  Tarawa  (LHA-1): 

•  Linear  defects  were  located  by  MT  inspection  in  original  fabrication  welds  and  re¬ 
pair  welds  alike,  including  welds  installed  as  recently  as  1  ‘*82. 

•  Predominant  linear  defect  was  the  transverse  weld  crack,  occurring  singly  and  in 
multipies  of  two  or  more;  crack  indications  ranged  from  1  16  inch  to  5  $  inch  in 
length. 

•  Numerous  linear  indications  were  also  found  across  the  thin  ligaments  between 
welds  in  the  A-to-D  loop  direction;  cracks  occurred  singly  and  in  multiples,  and 
ranged  in  length  from  1  16  inch  to  3  4  inch. 

«  Predominant  crack  orientation  was  transverse  to  the  longitudinal  axis  of  the 
header;  very  few  axially  oriented  indications  were  observed. 

•  Cracked  welds  in  the  inlet  outlet  header  were  concentrated  in  the  first  and  fourth 
passes. 

•  No  apparent  concentration  of  defects  in  the  welds  adjacent  to  the  division  plates 
(region  of  highest  thermal  gradient). 

•  The  inlet,  outlet  header  examined  contained  nearly  2-1  2  times  as  many  weld  joints 
with  linear  defects  as  did  the  intermediate  header  (32%  of  welds  in  no.  1  boiler 
inlet,  outlet  header  vs.  13%  of  welds  in  no.  2  boiler  intermediate  header). 

•  Destructive  analysis  indicated  the  maximum  crack  depth  to  be  approximately  5  32 
inch,  and  the  crack  propagation  mechanism  to  be  fatigue  at  elevated  temperature. 

Figure  l  on  page  3  and  Figure  2  on  page  4  are  photographs  taken  by  NAVSEA  show¬ 
ing  magnetic  particle  indications  of  the  cracking.  A  more  detailed  description  of  this 
problem  is  given  in  [Ref.  1:  pp.  10-15]. 
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II I.  DESCRIPTION  OF  ADINA 


ADINA  is  a  finite  element  program  for  automatic  dynamic  incremental  nonlinear 
analysis.  The  ADINA  system  is  composed  of  four  main  programs:  ADINA-IN, 
ADINA.  ADINA-PLOT  and  ADIN'AT.  The  1984  version  was  used  for  this  thesis.  The 
sequence  of  execution  of  these  programs  is  shown  in  Figure  3  on  page  6  which  is  taken 
from  [Ref.  2:  n.  1-2]. 

A.  GENERATION  OF  INPUT  DATA 

ADINA-IN  generates  the  input  data  for  ADINA  or  ADIN'AT.  ADINA  is  used  for 
displacement  and  stress  analysis  and  ADIN'AT  performs  the  temperature  analysis. 
ADINA-IN  can  be  run  in  an  interactive  or  batch  mode. 

Coordinates  can  be  input  in  a  local  or  global  coordinate  system  using  cartesian,  cy¬ 
lindrical  or  spherical  coordinates.  There  are  several  automatic  mesh  generation  methods 
available. 

ADINA-IN  has  the  ability  to  graphically  check  the  model.  All  or  part  of  the  mode1 
can  be  drawn  from  any  view  desired,  including  clement  and  nodal  numbering.  The 
LINES  parameter  described  in  [Ref.  2:  p.  5.25- 15J  gives  some  ability  to  specify  which 
lines  on  the  model  will  be  drawn.  This  allows  an  outline  of  the  model  to  be  drawn  and 
the  ability  to  check  for  missing  elements. 

The  procedure  for  utilizing  ADINA-IN  to  generate  input  data  for  ADINA  on  the 
VAX/VMS  operating  system  is  explained  in  Appendix  A.  The  procedure  for  utilizing 
ADINA-IN  to  generate  input  data  to  solve  the  thermally  induced  stress  problem  is  ex¬ 
plained  in  Appendix  B.  The  reverse  Cuthill-McKee  bandwidth  optimization  is  used  in 
generating  the  input  file  for  ADINA. 

13.  POST-PROCESSING 

ADINA-PLOT  performs  the  post-processing  of  the  output  data  from  ADINA.  It 
also  has  graphics  capability,  allowing  a  view  of  both  the  deformed  and  undeformed 
shapes.  The  deformed  scale  factor  can  be  specified.  ADINA-PLOT  can  also  do  vector 
plotting  of  the  principal  stresses. 

Several  listing  options  are  available.  Displacements  or  stresses  for  all  or  part  of  the 
model  can  be  printed.  Maximum  v.  lues  or  results  exceeding  a  specified  value  can  also 
be  listed. 


The  procedure  for  operating  ADINA-PLOT  on  the  VAX/VMS  operating  system  is 
explained  in  Appendix  C. 

C.  METHOD  OF  SOLUTION 

ADIXA  uses  isoparametric  displacement  based  finite  elements.  The  basic  assump¬ 
tions  from  [Ref.  3:  p.  2.3-8]  are: 


for  the  coordinates 


.=1  i=l  i=l 


for  the  displacements 


q 


llU[ 


;'-l 


v = Tm 


a 


i—  1 


(3.1) 


(3.2) 


where: 

•  /:,(r,s,t)  =  the  interpolation  function  corresponding  to  node  i 

•  r.s.t  =  the  isoparametric  coordinates 

•  q  =  the  number  of  element  nodes 

•  x,y„z,  =  nodal  point  coordinates 

•  u„ i’„w,  =  nodal  point  displacements 

ADINA  uses  the  principle  of  virtual  displacement  to  generate  the  equilibrium 
equations.  The  equilibrium  equation  as  listed  in  [Ref.  4:  p.  124J  is: 

KU=R  (3.3) 


where: 

•  K  =  the  clement  stiffness  matrix 

•  U  =  the  unknown  nodal  point  displacements 

»  R  =  the  load  vector  which  includes  body  forces,  surface  forces,  initial  stresses  and 
concentrated  loads. 


The  development  of  the  stiffness  matrix  for  the  isoparametric  clement  is  given  in  [Ref. 
4:  pp.  203-204]: 


A  =  B‘  C  B  dV 


where: 


•  B  =  the  strain-displacement  matrix 


•  C  =  the  elasticity  matrix 


In  terms  of  the  isoparametric  coordinates: 


d V  =  det(j)  dr  ds  dl 


where: 


det(J)  =  the  determinant  of  the  Jacobian  operator 


The  stresses  arc  calculated  using  the  strains  at  the  point  of  interest.  The  basic  relation¬ 
ships  used  are  given  in  [Ref.  4:  p.  123]: 


c  =  BU 


T  =  C  £ 


•  £  =  the  element  strain  matrix 


•  r  =  the  element  stress  matrix 


The  element  matrices  are  calculated  using  Gauss  numerical  integration.  A  (3x3x3)  in¬ 
tegration  order  is  utilized  as  specified  in  [Ref.  3  :  p.  2.3-16]. 


D.  AXISYMMETRIC  PROBLEM 


The  basic  assumptions  for  the  axisymmetric  element  are  listed  in  [Ref.  3:  p.  2.2-4]: 


£xx  —  ujy  (.v  is  i he  hoop  direction) 


rxy  —  rX2  =  0  (y  is  the  radial  direction) 


If  the  loadins  is  also  axisvmmetric,  a  2-dimcnsional  analysis  of  a  unit  radian  of  the 


structure  will  provide  the  complete  stress  and  strain  distributions.  This  greatly  reduces 
the  computational  time  compared  to  a  3-dimensional  analysis. 


E.  CONTACT  SURFACES 


ADINA  allows  the  specification  of  contact  surfaces  in  the  development  of  the 
model.  Contact  surfaces  are  defined  in  [Ref.  3:  p.  4.1-1]  as  surfaces  that  arc  initially  in 
contact,  or  are  anticipated  to  come  into  contact  during  the  response  solution. 

The  contact  surface  problem  is  a  nonlinear  analysis.  After  each  equilibrium  iter¬ 
ation,  the  segment  contact  force  is  evaluated  from  the  nodal  forces.  A  contactor  seg¬ 
ment  is  shown  in  Figure  4  on  page  10.  The  conditions  of  either  tension  release,  sliding 
or  sticking  are  determined  for  each  segment.  The  basic  equations  are  given  in  [Ref.  3: 
p.  4.2-3].  Tension  release  occurs  if: 

?.7^<0  (3.10) 


where: 

3  7  =  the  inward  normal  vector  to  segment  j  as  shown  in  Figure  4  on  page  10 

*  7;  =  the  normal  contact  force  for  segment  j 

Tension  release  means  the  surfaces  are  not  in  contact. 

If  the  surlaces  are  in  contact  the  condition  of  sliding  or  sticking  is  determined  by 
using  the  friction  law.  Sticking  occurs  if: 

(3.11) 


where: 

•  V,  =  the  tangential  contact  force  for  segment  j 

•  a  =  the  coefficient  of  friction 

Sliding  occurs  if: 


(3.12) 

The  nodal  forces  are  calculated  as  the  consistent  nodal  loads  which  correspond  to  the 
tangential  and  normal  forces  of  the  segment. 

F.  THERMAL  ANALYSIS 

ADI  NAT  uses  isoparametric  finite  element  discretization.  The  basic  assumptions 
from  [Ref.  5:  p.  11]  are: 
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(3.13) 


i  v  ft  't 

t  +  At9  =  Ya,  +  Ar„  ,  r*  =  TV.  i  =  YhAdi 


where: 


*  h.  =  the  element  interpolation  function 

*  N  =  the  number  of  nodal  points 

*  .v.j;,zr.  =  coordinates  of  the  nodal  point  i 

*  t  +  Ate ,  =  the  temperature  cf  nodal  point  i  at  time  t  +  At 
J  rv,  =  the  temperature  of  nodal  point  i  at  time  t 

»  Ad,  =  the  temperature  increment  at  nodal  point  i 


The  interpolation  and  substitution  of  these  equations  into  the  general  incremental  heat 
How  equilibrium  equations. yield  the  finite  element  equations. 


G.  VERIFICATION 


Part  of  the  installation  process  of  a  finite  element  code  is  the  verification  of  the  code. 
This  is  done  by  modeling  a  problem  for  which  the  results  are  known  and  comparing  the 
results. 


Three  verifications  are  included  in  this  thesis.  The  solution  of  a  cantilever  beam 


problem  is  compared  with  classical  beam  theory  in  Appendix  D.  The  thermal  stresses 
in  an  annular  disk  are  compared  with  classical  disk  theory  in  Appendix  E.  The  temper¬ 
ature  distribution  in  a  chimney  is  compared  with  a  verified  finite  difference  solution  in 
Appendix  F.  The  last  two  problems  were  chosen  because  of  their  similarity  to  the 
problem  addressed  in  this  thesis.  In  each  case  ADINA  vielded  excellent  results. 


IV.  2-DIMENSIONAL  SOLID  AXISYMMETRIC  FINITE  ELEMENT 

MODEL 


This  problem  was  originally  modeled  by  Lieutenant  Jon  W.  Kaufmann,  United 
States  Navy,  in  [Ref.  1].  The  initial  attempt  of  this  thesis  was  to  duplicate  his  model. 

A.  INITIAL  MODEL  DEVELOPMENT 

NAVSEA  provided  a  cross  section  drawing  of  the  superheater  header  as  illustrated 
in  Figure  5  on  page  13.  They  also  provided  a  detaii  of  the  superheater  header  tube  at¬ 
tachment  weld  joint  in  Figure  6  on  page  14.  The  same  geometry  of  the  weld  as  devel¬ 
oped  in  [Ref.  1:  p.  IS]  was  utilized  and  is  illustrated  in  Figure  7  on  page  15. 

A  nine-node  plane  axisymmetric  element  was  used.  This  element  is  similar  to  the 
one  used  in  [Ref.  1:  p.  29], 

A  coarse  and  fine  mesh  were  developed  in  order  to  do  a  study  on  convergence.  The 
coarse  mesh  has  64  elements  and  292  nodes.  The  model  is  shown  in  Figure  S  on  page 
16.  A  close-up  of  the  weld  is  shown  in  Figure  9  on  page  17,  listing  element  numbers. 
The  fine  mesh  has  224  elements  and  1095  nodes  and  is  illustrated  in  Figure  10  on  page 
13.  The  weld  area  is  shown  in  Figure  11  on  page  19. 

The  same  boundary  conditions  stated  in  [Ref.  1:  p.  30]  were  applied.  The  right  hand 
side  of  the  model  was  fixed  as  illustrated  in  Figure  12  on  page  20.  A  frictional  surface 
was  applied  between  the  tube  and  the  header.  The  nodes  on  the  boundary'  were  con¬ 
strained  together  in  the  y  direction  but  allowed  to  move  with  respect  to  each  other  in  the 
z  direction.  This  is  in  elTect  applying  a  contact  surface  with  the  coefficient  of  friction 
equal  to  zero. 

B.  LOAD  CONDITIONS 

The  same  load  conditions  were  applied  as  developed  in  [Ref.  1:  pp.  38-42]. 

1.  Load  condition  no.  1  -  Internal  pressure  load 

Lead  condition  no.  1  consisted  of  an  internal  pressure  of  700  psi  applied  at  all 
nodes  located  on  the  inside  of  the  tube  and  header.  This  corresponds  to  the  full  rating 
of  the  boiler.  This  loading  is  illustrated  in  Figure  13  on  page  21. 
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Figure  6.  Detail  of  the  attachment  weld  joint 
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2.  Load  condition  no.  2  -  Longitudinal  tube  load 

The  force  per  radian  is  152  Ibf.  The  coarse  model  uses  the  same  consistent 
loading  as  [Ref.  1:  p.  40].  The  load  was  applied  to  the  bottom  element  of  the  tube,  on 
nodes  S5,  68  and  S7  as  indicated  in  the  following  list: 

•  Fss  =  F  (1/6)  =  25.3  (Lb,) 

•  (2/3) -101.3  (Lb,) 

•  Fr  =  F  (1/6)  =  25.3  (Lb,) 

This  is  illustrated  in  Figure  14  on  page  23.  The  consistent  loading  for  the  fine  mesh  was 
applied  to  nodes  85,  240,  239,  238  and  87  as  follows: 

•  F3S  =  F  (1/12)  =  12.66666666  (Lb,) 

•  Fm  =  F  (1/3)  =  50.66666666  (Lb,) 

•  =  F  (1/6)  =  25.33333333  (Lb,) 

•  .Si  =  F  ( 1/3)  =  50.66666666  ( Lb,) 

»  />  =  F  (1/12)=  1 2.66666666  ( Lb.) 

The  consistent  loading  for  the  fine  mesh  is  illustrated  in  Figure  15  on  page  24. 

3.  Load  condition  no.  3  -  Thermal  load 

The  thermal  load  consists  of  maintaining  a  constant  temperature  of  350  degrees 
F  at  the  interior  nodes  of  the  tube  and  header.  A  temperature  gradient  of  100  degrees 
F  is  applied  to  the  right  hand  side  of  the  model.  This  load  condition  is  shown  in 
Figure  16  on  page  25. 

C.  MATERIAL  PROPERTIES 

The  same  material  properties  were  used  as  listed  in  [Ref.  1:  p.  3SJ.  Young's  Modulus 
(E)  equals  29,600,000  psi,  the  coellicient  of  thermal  expansion  (a)  equals  .0000065  in,  in, 
°F  and  Poisson's  ratio  (v)  equals  .3.  The  material  of  the  tube  and  header  is  2.25% 
chromium.  1%  molybdenum  steel  (ASTM  grade  22).  The  weld  was  assumed  to  have  the 
same  properties  as  the  tube  and  header. 

D.  COMPARISON  OF  RESULTS 

The  original  and  deformed  shapes  for  load  conditions  1,  2  and  3  for  the  coarse  mesh 
are  shown  in  Figure  17  on  page  26,  Figure  IS  on  page  27,  and  Figure  19  on  page  28. 
Tiie  results  for  the  fine  mesh  are  shown  in  Figure  20  on  page  29.  Figure  21  on  page 
30,  and  Figure  22  on  page  31.  The  maximum  stresses  in  the  weld  for  the  coarse  and  fine 
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Figure  18.  Original  and  deformed  shape  of  coarse  mesh  for  load  condition  no. 
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Figure  19.  Original  and  deformed  shape  of  coarse  mesh  for  load  condition  no.  3 
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Figure  20.  Original  and  deformed  shape  of  fine  mesh  for  load  condition  no.  I 
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Figure  21.  Original  and  deformed  shape  of  fine  mesh  for  load  condition 
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Figure  22.  Original  and  deformed  shape  of  fine  mesh  for  load  condition 


mesh  are  listed  in  Tabie  1  and  Table  2.  For  the  axisvmmetric  problem,  is  the  hoop 
stress,  oJf  is  the  radial  stress,  and  oa  is  the  longitudinal  stress. 


Table  I.  MAXIMUM  STRESSES  FOR  THE  COARSE  MESH 


Stresses 

Load  Condition 

No.  1 

Load  Condition 

No.  2 

Load  Condition 

No.  3 

o„  (psi) 

-.636767E  +  03 

-.509766E  +  03 

-.135183E  +  06 

<7,v  ( psi) 

-.66 1331  £  +  03 

-.370531E  +  03 

-.466026E  +  05 

a.t  (psi) 

-.742S39E  +  03 

-.1 16857E-‘-03 

-.55943“E  — 04 

a,,  tpsi) 

.  1 209  3  3  E  —  0  3 

.332004E  +  03 

.  1 34600E  +  05 

Table  2.  MAXIMUM  STRESSES  FOR  THE  FINE  MESH 


Stresses 

Load  Condition 

No.  1 

Load  Condition 

No.  2 

Load  Condition 

No.  3 

ex„  i  psi) 

-.649740E  +  03 

-.534142E  +  03 

-.136960E  +  06 

! psi) 

.5035S7E-03 

-.555045E-05 

t_  (psi) 

145094E  —  03 

-.9S0113E  +  04 

n  .  <  psi  ) 

.  I3647SE  —  03 

.37294SE-03 

.  I93084E  — 05 

The  magnitude  and  direction  of  the  principal  stresses  for  the  coarse  and  fine  mesh, 
in  the  area  of  the  weld,  are  listed  in  Table  3  and  Table  4.  A  vector  plot  of  the  principal 
stresses  of  the  coarse  and  fine  mesh  for  load  condition  no.  3  is  shown  in  Figure  23  on 
page  34  and  Figure  24  on  page  35. 


Tahle  3.  PRINCIPAL  STRESSES  FOR  THE  COARSE  MESH 


Element 

Pi  (psi) 

P2  (psi) 

Angle  (°F) 

12 

-.1594E-03 

-.3422E  +  05 

74.57 

13 

0. 1352E  —  04 

l820E-f  05 

85.39 

14 

0.2S02E  +  03 

-,14I3E-r05 

94.46 

15 

0.9 1 53  B  —  03 

-.3779E  +  05 

86.39 

16 

0.3717E-04 

-.2S07E  +  05 

86.38 

17 

0.2304E  +  04 

86.18 

18 

0.1716E^04 

-.3S08E  +  05 

S8.03 

19 

0.3323E  +  04 

-.3259E  +  05 

87.72 

The  maximum  hoop  stress  listed  in  [Ref.  1:  pp.  39-40]  for  load  condition  no.  1  is  .69 
(ksi)  and  the  maximum  hoop  stress  for  load  condition  no.  2  is  .96  (ksi).  Those  results 
compare  favorably  with  the  results  shown  here  from  ADINA.  [Ref.  1:  p.  42]  indicates 
the  maximum  hoop  stress  for  load  condition  no.  3  is  19.0  (ksi).  This  value  is  much 
smaller  than  the  result  of  137  (ksi)  obtained  with  ADINA.  In  considering  the  verifica¬ 
tion  of  ADINA  as  explained  in  chapter  three,  it  is  believed  that  the  results  obtained  here 
are  correct.  The  results  confirm  that  the  thermal  stresses  are  the  predominant  stresses. 
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Figure  23.  Vector  plot  of  the  principal  stresses  of  the  coarse  mesh  for  load  condition 
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Figure  24.  Vector  plot  of  the  principal  stresses  of  the  fine  mesh  for  load  condition 


The  accuracy  of  the  results  of  the  problem  that  has  been  modeled  are  very  good. 
There  is  very  little  difference  between  the  results  of  the  coarse  and  fine  mesh  which 
means  convergence  has  been  achieved.  There  is  no  significant  difference  in  the  calcu¬ 
lated  stresses  at  the  boundaries  of  adjacent  elements. 

The  next  step  was  to  use  the  contact  surface  option  of  ADINA  for  the  boundary 
between  the  tube  and  the  header.  This  allows  an  indication  of  the  coefficient  of  friction 
and  a  determination  of  its  effect  on  the  solution.  The  maximum  stresses  are  shown  in 
Table  5  and  Table  6  for  the  coarse  and  fine  mesh  for  a  coefficient  of  friction  of  zero. 
These  results  compare  very  well  with  the  previous  results. 


Table  5.  MAXIMUM  STRESSES  FOR  THE  COARSE  MESH,  CON¬ 
TACT  SURFACE  OPTION 


|  Stresses 

Load  Condition 

No.  1 

Load  Condition 

No.  2 

Load  Condition 

No.  3 

|  (psi) 

-,64o923E  +  03 

-.514228E  +  03 

-.1352651:4-06 

f7 ,  ( psi  ) 

-.05  ;  6141;  +  03 

.362.5591;  +  03 

-.46.523  i  F:  4-  05 

1  n~.  1  psi) 

-.747963  E  +  03 

-.1 12506E  +  03 

-.5586 13  E  4-  04 

!  ff„(p si) 

.1343S2E  +  03 

.3285.54E  +  03 

.134376E4-05 

Table  6.  MAXIMUM  STRESSES  FOR  THE  FINE  MESH,  CONTACT 
SURFACE  OPTION 


Stresses 

Load  Condition 

No.  1 

Load  Condition 

No.  2 

Load  Condition 

No.  3 

axx  (psi) 

-.646745  E +  03 

-.537235E  +  03 

-.  1 37053  E  + 06 

(Psi) 

-.657250E  +  03 

.495395E  +  03 

-.553930E  +  05 

n.,  (psi) 

-.75704  /  E  +  03 

-.  1 39909  E  + 03 

-.97S366E  +  04 

n,.  (psi) 

.139310E  +  03 

.  366954  E  +  03 

.  192709  E  + 05 

E.  FURTHER  MODEL  DEVELOPMENT 

A  detailed  study  of  the  problem  and  the  results  obtained  in  the  previous  section  in¬ 
dicated  some  changes  were  needed  in  the  model. 

It  was  discovered  that  the  model  as  developed  was  not  a  repeating  section.  The 
distance  between  the  center  of  the  tubes  is  two  inches.  The  right  hand  side  of  the  model 
is  two  inches  from  the  center  of  the  tube.  A  new  model  was  developed  with  the  right 
hand  side  being  one  inch  from  the  center  of  the  tube.  The  coarse  mesh  consists  of  43 


elements  and  211  nodes.  The  fine  mesh  has  172  elements  and  763  nodes.  The  models 
are  shown  in  Figure  25  on  page  38  and  Figure  26  on  page  39. 

Figure  19  on  page  2S  and  Figure  22  on  page  31  prompted  a  change  in  boundary 
conditions.  The  deformed  shapes  showed  that  the  header  was  not  allowed  to  expand 
uniformly  in  the  axial  direction.  The  bottom  node  of  the  right  hand  side  of  the  header 
was  fixed  and  the  rest  of  the  nodes  on  the  right  hand  side  were  allowed  to  roll  in  the  z 
direction.  These  boundary  conditions  are  shown  in  Figure  27  on  page  40. 

The  original  and  deformed  shapes  for  load  conditions  1,  2  and  3  for  the  coarse  mesh 
are  shown  in  Figure  2S  on  page  41.  Figure  29  on  page  42  and  Figure  30  on  page  43. 
Figure  30  on  page  43  shows  that  the  header  is  now  expanding  uniformly  in  the  axial 
direction.  The  maximum  stresses  are  indicated  in  Table  7  and  Table  S  for  the  coarse 
and  fine  mesh.  The  ADINA-IN  input  file  for  these  models  is  listed  in  Appendix  G. 


Table  7.  MAXIMUM  STRESSES  FOR  TIIE  COARSE  MESH,  RE¬ 
PEATING  SECTION 


Stresses 

Load  Condition 

No.  1 

Load  Conditions 

1  and  2 

Lead  Condition 

No.  3 

( psi ) 

-.42038  IE +  03 

-.566508E  +  03 

-.997953  E  + 05 

<7,v  (psi) 

-.597283E  +  03 

.792176E  +  03 

-.373304E  +  05 

o.._  (psi) 

-.825134E  +  03 

-.1 12S52E  +  04 

-.454326E  +  04 

( psi) 

■iMIriUMIf 

.838120E  +  03 

.104692E  +  05 

Table  8.  MAXIMUM  STRESSES  FOR  THE  FINE  MESH,  REPEATING 
SECTION 


Stresses 

Load  Condition 

No.  1 

Load  Conditions 

1  and  2 

Load  Condition 

No.  3 

(psi) 

-.473741 E  + 03 

-.595979E  +  03 

-.101 403  E  +  06 

<7,t.  ( psi) 

-.84O852E  +  03 

.  1 303  59  E  +  04 

-.447373E  +  05 

n..  ( psi ) 

-.922094E  +  03 

-.126150E  +  04 

-.796674E  +  04 

n,.  (psi) 

.3  13721 E  +  03 

.  102645E  +  04 

.1516S4E  +  05 

F.  THERMAL  GRADIENTS. 

The  previous  change  in  the  model  made  the  thermal  gradient  in  the  area  of  the  weld 
very  fiat.  It  was  desired  to  determine  what  effect  different  thermal  gradients  on  the  right 
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Figure  26.  Fine  mesh,  repeating  section 
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Figure  29.  Original  and  deformed  shape  of  coarse  mesh  for  load  condition  no.  2  , 
repeating  section 
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hand  side  of  the  model  would  have  on  the  results.  The  maximum  stresses  for  a  thermal 
gradient  of  100,  10,  2  and  0  degrees  F  are  shown  in  Table  9. 


Table  9.  MAXIMUM  STRESSES  FOR  DIFFERENT  THERMAL  GRADI¬ 
ENTS 


Thermal 

Gradient 

(psi) 

(Ps0 

(psi) 

oi2  ( psi) 

P'DT 

-.9979531; +  05 

UM 

-.454326E  +  04 

.  104692E  +  U5 

lnaF 

-.99S90SE  +  05 

-.373122E-05 

-.45401 7E  + 04 

.  10435SE  +  05 

2°F 

-.99S993E  +  05 

■B 

-.453990E  +  04 

.10432SE  +  05 

|  0 °F 

-.99901 5E  + 05 

-.373102E  +  05 

-.453983E  +  04 

.  10432 1 E  +  05 

Tire  results  were  surprising.  The  size  of  the  gradient  had  very  little  effect  on  the 
magnitude  of  the  stresses.  The  stresses  remained  very  high  even  for  a  zero  thermal  gra¬ 
dient.  A  study  of  these  results  indicated  the  need  for  different  boundary  conditions. 
Restricting  the  right  hand  side  of  the  model  to  move  in  the  z  direction  only  was  not  al¬ 
lowing  the  header  to  expand  outwardly  in  the  radial  direction.  The  header  was  being 
forced  to  expand  into  the  center  of  the  tube  which  was  creating  the  high  stresses  in  the 
area  of  the  weld. 

The  boundary  condition  chosen  is  shown  in  Figure  31  on  page  45.  The  bottom  left 
hand  node  of  the  header  at  the  boundary  of  the  tube  and  the  header  was  fixed.  The  rest 
of  the  bottom  nodes  of  the  header  were  allowed  to  roll  in  the  y  direction.  The  contact 
surface  option  of  ADINA  could  not  be  used  for  this  solution.  The  nodes  on  the 
boundary  were  constrained  together  in  the  y  direction  and  allowed  to  move  with  respect 
to  one  another  in  the  z  direction,  modeling  a  contact  surface  with  zero  coefficient  of 
friction.  This  was  with  the  exception  of  the  bottom  left  hand  node  of  the  header. 

An  iterative  procedure  was  then  followed  to  determine  the  solution  of  the  problem. 
After  each  run  of  the  problem,  the  bottom  left  hand  node  of  the  header  was  given  a 
permanent  displacement  to  the  right  equal  to  the  displacement  of  the  corresponding 
node  of  the  tube  at  the  boundary.  This  was  continued  until  no  overlap  occurred.  The 
original  and  deformed  shapes  of  the  final  iteration  for  load  condition  no.  3  are  shown  in 
Figure  32  on  page  46.  As  can  be  seen,  the  model  is  now  expanding  uniformly  in  the 
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Figure  31.  Bottom  rolling  boundary  condition 


REPEAT  HEADER 

ORIGINAL  i _  . 

DEFORMED  t _ 

TIME  2.000 


D.34  9 
D.  00369 


-  y 


f- — -j-v; 


■  i  !  >  :  i  ! 
r - H  !  | 


i  n  1 5 


i  i  I  I!  j  !  I 

:  I'tlU 


""77  i 

mi: 


i  i 

*  I  ! 


Figure  32.  Original  and  deformed  shape  of  the  fine  mesh  for  load  condition  no.  3. 
expanding  radially  and  axiallv 
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axial  direction  and  outwardly  in  the  radial  direction.  The  maximum  stresses  in  the  weld 
area  of  the  coarse  -nosh  for  load  condition  no.  3  are  shown  in  Table  10. 


Table  10.  MAXIMUM  STRESSES  FOR 
THE  COARSE  MESH,  EXPAND- 


ING  RADIALLY 

Stresses 

Load  Condition  No.  3 

( psi) 

.460123E  +  03 

a  .  ( psi > 

-.932915E  +  02 

:  !  psi ) 

.1022271-:  4-  03 

*  psi » 

-.584215E  +  02 

These  results  are  much  lower  than  those  obtained  for  the  previous  model  and  it  is 
believed  that  these  results  are  closer  to  the  actual  stresses.  However,  it  is  also  believed 
that  this  model  is  not  a  true  repeating  section  in  that  the  expansion  of  the  neighboring 
tubes  and  the  rest  of  the  header  will  have  an  effect  on  this  section,  probably  raising  the 
stresses  in  the  weld  area.  The  need  for  a  3-dimensional  model  was  indicated. 

A  3-dimensionai  model  was  also  desired  in  order  to  apply  non-symmctric  loading  to 
the  model.  The  major  difference  between  this  header  and  other  headers  with  a  similar 
cross  section  is  that  this  header  is  not  sectionalized.  It  has  one  inlet; outlet  header  and 
one  intermediate  header.  This  results  in  a  larger  temperature  gradient  from  the  inlet  to 
the  outlet  side.  A  3-dimensional  model  would  allow  this  type  of  non-symmetric  loading 
to  be  studied. 


V.  3-DIMENSIONAL  SOLID  FINITE  ELEMENT  MODEL 

A.  MODEL  DEVELOPMENT 

The  repeating  section  chosen  to  be  modeled  was  one  tube  row.  Only  half  of  the  row 
was  modeled  due  to  symmetry.  The  welds  were  assumed  to  have  the  same  material 
properties  as  the  tube  and  header  and  no  attempt  was  made  to  substructure  the  weld 
area. 

The  element  chosen  to  model  the  structure  was  the  20-node  isoparametric  solid  ele¬ 
ment.  This  element  is  shown  in  Figure  33  on  page  49. 

Figure  34  on  page  50  shows  the  modeling  of  one  tube.  Each  tube  consists  of  sixty 
elements.  Figure  35  on  page  51  shows  how  the  heudei  was  connected  to  the  tube.  A 
contact  surface  was  placed  between  the  tube  and  the  header.  The  entire  model  is  shown 
in  Figure  36  on  page  52.  Figure  37  on  page  53  is  another  view  of  the  model  showing 
the  location  of  the  origin  of  the  axis.  The  model  consists  of  2SS  elements  and  1957 
nodes.  Figure  3S  on  page  54.  Figure  39  on  page  55  and  Figure  40  on  page  56  show  a 
view  of  the  model  looking  from  the  x,  y  and  z  directions. 

The  hidden  line  computer  code  developed  by  David  R.  Iledgley.  Jr.  of  the  Ames 
Research  Center  was  utilized  to  check  the  model.  Different  views  are  shown  in 
Figure  41  on  page  57  and  Figure  42  on  page  58. 

The  following  boundary  conditions  were  imposed.  Every  node  in  the  y  =  0  plane 
was  constrained  to  roll  in  that  plane.  In  addition  the  node  at  the  origin  was  fixed  and 
the  nodes  above  the  origin  along  the  positive  x  axis  were  constrained  to  roll  only  in  the 
x  direction. 

B.  INTERNAL  PRESSURE  LOAD 

The  internal  pressure  load  of  700  psi  was  applied  to  the  inside  nodes  of  the  tube  and 
header.  It  was  then  discovered  that  ADINA  does  not  optimize  the  band-width  of  con¬ 
tact  surfaces.  The  bandwidth  of  this  model  is  approximately  2700  and  the  problem  was 
too  large  to  be  solved  on  the  VAX.  VMS  operating  system.  The  ADINA-IN  input  file 
for  the  contact  surfaces  is  listed  in  Appendix  II  for  reference. 

A  contact  surface  was  then  modeled  by  constraining  the  corresponding  nodes  of  the 
tube  and  header  to  move  together  in  the  y  and  z  direction,  but  allowing  ti.'tn  to  move 
relative  to  one  another  in  the  x  direction.  This  assumes  a  coefficient  of  frictic.i  of  zero 
and  reduced  the  bandwidth  to  approximately  1000. 
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Figure  38.  View  of  3-dimensional  model 


It  then  became  necessary  to  enforce  symmetry'  on  either  side  of  the  z  =  0  plane. 
Every  node  in  this  plane  was  constrained  to  roll  in  the  plane.  It  is  understood  that  this 
enforced  symmetry  is  more  costly  than  modeling  only  half  of  the  structure.  The  entire 
section  was  modeled  to  allow  a  study  of  unsymmetric  temperature  loading.  The 
ADINA-IN  input  file  is  listed  in  Appendix  I. 

Figure  43  on  page  60  shows  the  original  and  deformed  shapes  for  this  loading  con¬ 
dition.  Table  1 1  lists  the  maximum  stresses  in  the  area  of  the  weld. 


Table  11.  MAXIMUM  STRESSES 


3-DIMENSIONAL  MODEL 


|  Stresses 

Internal  pressure  load 

!  ipsi) 

-.9S4072E  +  03 
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Figure  43. 
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VI.  CONCLUSIONS 


A.  DISCUSSION  OF  RESULTS 

The  maximum  stresses  obtained  for  the  axisymmetric  model  were  between  137  (ksi) 
and  .46  (ksi),  dependent  upon  the  boundary'  conditions  imposed.  The  axisymmetric 
model  does  not  adequately  represent  the  physical  boundary  conditions  and  the  actual 
stress  is  somewhere  between  these  two  figures. 

A  3-dimensional  model  has  been  developed  which  allows  an  accurate  representation 
of  the  boundary  conditions.  The  maximum  stress  for  the  internal  pressure  load  in  the 
area  of  the  weld  was  1.S5  (ksi).  This  model  also  allows  non-symmctric  loading  to  be 
applied. 

3.  OPPORTUNITIES  FOR  FURTHER  RESEARCH 

The  next  step  in  solving  this  problem  is  to  obtain  results  for  the  thermal  loading. 
The  initial  loading  should  be  static  but  the  transient  loading  should  also  be  studied.  In 
studying  the  transient  problem,  a  consideration  should  be  made  of  the  temperature  de¬ 
pendence  of  the  material  properties  of  the  structure.  Another  area  for  research  would 
be  the  consideration  of  creep. 

The  major  attempt  of  this  report  was  to  develop  a  3-dimensional  model  which  could 
be  used  to  investigate  the  failures  that  have  occurred  in  the  superheater  header  attach¬ 
ment  welds.  This  model  can  also  be  used  to  investigate  any  needed  design  changes.  It 
is  hoped  that  this  model  will  be  used  for  further  investigation  of  this  problem. 


APPENDIX  A.  ADINA-IN  PROCEDURES 


These  procedures  are  for  the  VAX,  VMS  operating  system  utilizing  the  plot-10 
graphics  library. 

To  operate  ADINA-IN  in  the  batch  mode,  a  data  file  which  defines  the  problem 
must  be  created.  Assign  this  data  file  to  a  logical  unit  number  and  run  ADINA-IN. 

Now  create  the  database.  All  of  the  input  to  ADINA-IN  is  stored  on  this  database. 
The  Command  is  DATAB  CREATE.  Next  read  the  data  file  by  using  the  command: 


READ  (logical  unit  number). 

The  model  can  now  be  plotted  using  the  frame,  mesh  and  plot  commands  as  ex¬ 
plained  in  [Ref.  2:  pp.  5.25-1  to  5.25-20J.  The  CONTROL  command  must  be  set  to  the 
correct  parameters  to  use  graphics.  Tor  the  VAX  VMS  operating  system  utilizing  the 
plot-10  graphics  library,  set  buudratc  =  4S00,  system=2,  device  =  2  and  option  =  1024. 


Chances  can  be  made  to  the  model  bv  cxitins  and  chancing  the  data  file  or  bv  tvpinu  in 


changes  interactively.  Ail  changes  will  be  saved  on  the  database. 

When  the  problem  is  properly  defined,  create  the  ADINA  input  file  using  the  com¬ 
mand  ADINA.  Exit  ADINA-IN  by  typing  in  END. 

ADINA-IN  will  have  created  two  files.  FOR001.DAT  is  the  database.  This  data¬ 


base  can  best  be  saved  by  renaming  it  to  a  unique  file  name.  To  make  changes  assign 
the  database  to  logical  unit  FOROOl  and  run  ADINA-IN.  Use  the  command  DATA¬ 


BASE  OPEN  and  then  make  any  necessary  changes. 

The  second  file  created  by  ADINA-IN  is  FOR002.DAT  and  this  is  the  input  file  for 
ADINA.  This  file  should  be  renamed  to  a  unique  filename  and  assigned  to  logical  unit 
number  FOR005.  Assign  an  output  file  to  logical  unit  number  FOR006.  Ensure  there 
are  no  FOR"'."'  files  in  the  directory  as  this  can  interfere  with  the  operation  of  ADINA. 
Now  run  ADINA.  The  requested  results  will  be  printed  in  the  assigned  output  file. 


APPENDIX  B.  ADINA-IN  PROCEDURES  FOR  THE  THERMAL  STRESS 

PROBLEM 


ADINA-IN  is  also  used  for  creating  the  input  file  for  ADINAT.  This  particular 
version  of  ADINA-IN  will  not  create  the  entire  file  which  necessitates  some  minor  ed¬ 
iting  of  the  ADINAT  input  file  before  ADINAT  can  be  run. 

The  temperature  loading  in  this  thesis  was  done  by  using  the  LOADS  TEMPER¬ 
ATURE  command  in  the  ADINA-IN  input  file.  After  the  problem  is  defined  in 
ADINA-IN.  use  the  ADINAT  command  to  generate  the  ADINAT  input  file. 

ADINA-IN  again  creates  two  files.  FOR001.DAT  is  the  database  and 
FOROG2.DAT  is  the  ADINAT  input  file.  Rename  these  files. 

The  ADINAT  input  file  must  now  be  edited.  Search  for  the  line  which  has  the  let¬ 
ters  'KKKKKK'.  This  line  must  be  deleted  and  the  specific  heat  entered.  The  line  just 
above  this  line  must  also  be  deleted  and  the  thermal  conductivity  entered.  Exit  the 
ADINAT  input  file. 

Assign  the  ADINAT  input  file  to  logical  unit  number  FOROU5  and  assign  an  output 
file  to  logical  unit  number  FGR006.  Run  ADINAT.  ADINAT  creates  several  files, 
most  of  which  are  used  for  the  operation  of  the  code.  The  final  nodal  point  temperatures 
and  heat  fluxes  are  printed  in  the  output  file.  The  nodal  temperatures  are  also  stored  in 
FOR056.DAT.  This  file  must  be  renamed  to  a  unique  filename. 

ADINA-IN  must  now  be  run  again  to  create  the  input  file  for  ADINA  to  solve  for 
the  thermal  stresses.  Assign  the  nodal  temperature  file  to  logical  unit  number  FOR056. 

The  ADINA-IN  input  file  must  be  changed.  Delete  the  LOADS  TEMPERATURE 
command  and  indicate  on  the  MASTER  command  line  that  ITP56=  1.  Run 
ADINA-IN  and  read  the  changed  input  file.  Ensure  all  FOR*.*  files  are  deleted  before 
doing  this.  Then  use  the  ADINA  command  to  create  the  ADINA  input  file.  ADINA 
can  now  be  run  as  explained  in  the  last  paragraph  of  Appendix  A. 


APPENDIX  C.  ADINA-PLOT  PROCEDURES 


ADINA  created  several  files  besides  the  assigned  output  file.  Most  of  these  files  are 
used  for  the  operation  of  the  code.  File  FOR060.DAT  is  the  porthole  file  for 
ADINA-PLOT. 

Before  running  ADINA-PLOT,  ensure  all  FOR001.DAT  files  are  deleted  and  that 
logical  unit  number  FORGO  1  is  deassigned.  Run  ADINA-PLOT.  First  create  the  da¬ 
tabase  using  the  DATABASE  CREATE  command.  All  of  the  information  in  the 
porthole  file  is  now  loaded. 

Ensure  the  control  information  is  correct.  The  model  can  now  be  plotted  and  the 
results  listed  as  explained  in  [Ref.  6). 

Exit  ADINA-PLOT  with  the  command  END.  ADINA-PLOT  will  have  created  a 
database  with  filename  FOR.OOl.DAT.  Rename  this  file.  To  look  at  the  results  again, 
assign  the  database  to  logical  unit  number  FORGO  1  and  use  the  DATABASE  OPEN 
command  when  running  ADINA-PLOT. 


APPENDIX  D.  VERIFICATION  OF  ADINA  -  CANTILEVER  BEAM 

PROBLEM 


The  cantilever  beam  problem  modeled  is  shown  in  Figure  44  on  page  66.  A  moment 
of  360  k-in  is  applied  to  the  end  of  a  cantilever  beam  which  is  48  inches  long  and  12 
inches  wide. 

The  classical  equations  as  listed  in  [Ref.  7:  pp.  215  and  737]  are: 


a  =  My!  I 

(.4.1) 

b  -  M0L2I2EI 

(A.2) 

The  finite  element  model  showing  both  the  original  and  deformed  shapes  is  shown  in 
Figure  45  on  page  67.  A  nine  node  plane  stress  element  was  used.  Table  12  shows  a 
comparison  of  the  results  for  several  nodes.  As  can  be  seen,  ADINA  yields  exact  results 
for  this  problem.  The  following  is  the  ADINA- IN  input  file: 


Table  12.  COMPARISON  OF  RESULTS 


Comparison 

Beam  Theory 

ADINA 

displacement  of  node  10  (inch) 

.096 

.096 

stress  at  node  1  (psi) 

15 

15 

stress  at  node  12  (psi) 

-15 

-15 

DATAB  CREATE 

HEAD  '  CANTILEVER  BEAM' 

A 

MASTER  REACTIONS = YES  ID0F=100111 
PRINTOUT  MAX 
ANALYSIS  TYPE=STATIC 

A 

COORDINATES  /  ENTRIES  NODE  X  Y  Z 

1  0  0  0 

5  0  48  0 

15  0  48  12 

11  0  0  12 

A 


I 

i  65 

\ 


> 


Ml 


MATERIAL  1  ELASTIC  E=3E4  NU=0. 0  D=1 

A 

EGROUP  1  PLANE  SUBTYPE = STRESS 2  RESULTS = STRESSES  INT=6 
STRES STABLE  N=1  Pl=l  P2  =  2  P3=3  P4=4  P5=5  P6=6  P7=7  P3=8  P9=9 
EDATA  /  ENTRIES  EL  THICK 

1  1 
2  1 

GSURFACE  1  5  15  11  ELI =2  EL2=1  NODES =9 

A 

BOUNDARIES  111111  TYPE=NODES  /  1  6  11 

A 

LOADS  CONCENTRATED 
5  2  30 
15  2  -30 

A 

CONTROL  B=4800  SYSTEM=2  DEVICE=2  0PTI0N=1024 
CONTROL  PLOTE=W  PLOTS =N 


APPENDIX  E.  VERIFICATION  OF  ADINA  -  ANNULAR  DISK 

PROBLEM 


The  classical  equations  for  the  thermal  stresses  are  listed  in  [Ref.  8  :  p.  259]: 


or  =  a  £[  —  1 


Trdr  4-  (r2  —  a2)lr'(b2 


(EA) 


=  aET  —T+  1  jr 


Trdr  +  (r2  +  a~)/r2(b2  —  a2)  Trdr ] 

*  n 


(• E.2 ) 


where: 

•  a  =  inner  radius,  .63  inch 

•  b  =  outer  radius,  3.0  inch 

•  T  =  temperature  distribution,  399.569212  -  86. 3436r  +  12.16239S67H 

•  E  =  29600000  psi 

•  a  ~  .0000065  in/in;  °F 

The  finite  element  model  showing  the  original  and  deformed  shapes  is  shown  in 
Figure  46  on  page  70.  A  nine  node  plane  axisymmetric  element  was  used.  Table  13 
shows  a  comparison  of  the  results.  The  comparison  is  very  good.  Most  of  the  difference 
can  be  attributed  to  the  fact  that  the  results  for  the  disk  theory  were  computed  at  the 
nodes  where  as  the  results  for  ADINA  are  calculated  at  the  Gaussian  integration  points. 
The  following  is  the  ADINA- IN'  input  file: 


Table  13.  COMPARISON  OF  RESULTS 


r  (inch) 

Stresses  from  Disk  Theory 
(psi) 

Stresses  from  ADINA 
(psi) 

.63 

o„  —  -.139E  +  05 

Off  =  -.  120961 E +  05 

1.815 

=  -.2577E  +  04 

<7,  =  -.253146E  +  04 

=  .  17I5E  +  04 

o„  —  .  169057E  +  04 

_ 

0^  =  .534E  +  04 

Off  =  .51 1 644  E +  04 

DATAB  CREATE 
HEAD  'DISC' 

CONTROL  B=4800  SYSTEM=2  DEVICE=2 
CONTROL  PLOTE=WAIT  PLOTS=NO 


OPTION=1024 


COORDINATES  /  ENTRIES  NODE  X 

1  0 

2  0 

3  0 

4  0 

A 


Y  Z 
.63  0 

.63  .3 

3.0  0 

3.0  .3 


MATERIAL  1  THERMO-ELASTIC 
-2000  29600000  .3  .0000065 

2000  29600000  .3  .0000065 


EGROUP  1  PLANE  SU=AXISYMM 
GSURFACE  4213  ELI =4  EL2=1  N0=9 

A 


MASTER  REACTIONS = YES  IDOF*100111 
PRINTOUT  MAX 

ANALYSIS  TYPE = STATIC 
* 


BOUNDARIES  101111  TYPE  =N0DES/12  STEP  1  TO  20 
A 

LOADS  TEMP  TR=0 

1  350 

2  350 
20  350 

3  250 

4  250 
12  250 
* 
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APPENDIX  F.  VERIFICATION  OF  ADINA  -  CHIMNEY  PROBLEM 


The  chimney  problem  modeled  is  shown  in  Figure  47  on  page  73.  The  inside  surface 
is  at  a  uniform  temperature  of  650  degrees  C  and  the  outside  surface  is  at  a  temperature 
of  150  degrees  C. 

The  finite  element  model  is  shown  in  Figure  48  on  page  74.  A  nine  node  plane 
stress  element  was  used. 

A  finite  difference  temperature  distribution  of  this  problem  was  obtained  from  Pro¬ 
fessor  G.  Cantin  of  the  Naval  Postgraduate  School.  The  results  had  converged  to  the 
exact  solution.  Table  14  shows  a  comparison  of  the  results  for  several  nodes.  The 
comparison  is  very  good.  The  following  is  the  ADINA-IN  input  file: 

Table  14.  RESULTS  FOR  CHIMNEY  PROBLEM 


Node 


Temperatures  from  Finite 
Difference  t°C) _ 

3.6933026E  +  02 

3.466 1577E  +  02 _ 

2.5149076E  +  02 _ 

3.4S91941E  +  02 _ 

3.8732346E-i*02 


Temperatures  form 
ADINA  ! °C) 

3.6S497E  -r  02 

3.39891  E  +  02 

2.52863E  +  02 

3.42042E  ■+•  02 

4.01579E  +  02 


DATAB  CREATE 
HEAD  'CHIMNEY' 

CONTROL  B=4800  SYSTEM=2  DEVICE=2  OPTION=1024 
CONTROL  ?L0TE=WAIT  PL0TS=NO 

h 

COORDINATES /ENTRIES  NODE  X  Y  Z 


2.0  1.0 


MATERIAL  1  THERMO -ELASTIC 
-2000  29600  .3  .0000065 
2000  29600  .3  .0000065 

* 


CH1MNE 

ORIGIN 


EGROUP  1  PLANE  SU=STRESS2 
EDATA 


1  1 
2  1 
3  1 

GSURFACE  1234  EL1=1  EL2=1  N0=9 


GSURFACE  4567  ELI =2  EL2=1  NO =9 

A 


MASTER  REACTIONS = YES  IDOF=100111 
PRINTOUT  MAX 
ANALYSIS  TYPE = STATIC 

A 

LOADS  TEMPERATURE 


APPENDIX  G.  ADINA-IN  INPUT  FILE  FOR  THE  AXISYMMETRIC 

PROBLEM 


OATAfi  CREATE 

* INPUT  FOR  THE  COARSE  MESH,  2-DIMENSIONAL  AXISYMMETRIC  MODEL 
HEAD  'REPEAT  HEADER' 

CONTROL  B-4800  SYSTEM-2  DEVICE-2  OPTION-1024 
CONTROL  PLOTE-HAIT  PLOTS-NO 
^COORDINATES  ARE  DEFINED 

COORDINATES  /  ENTRIES  NODE  X  Y  2 


1 

0 

.63 

3.0 

3 

0 

.693 

3.0 

37 

0 

.63 

2.844 

39 

0 

.75 

2.844 

60 

0 

.63 

2.782 

61 

0 

.75 

2.782 

70 

0 

.63 

1.0 

71 

0 

.75 

1.0 

72 

0 

.75001 

1.0 

85 

0 

.63 

0.0 

87 

0 

.75 

0.0 

100 

0 

.938 

3.0 

200 

0 

.82 

2.837 

152 

0 

.938 

2.97333333 

204 

0 

.9322366 

2.93059481 

211 

0 

.901 

2.88 

210 

0 

.84513864 

2.84814711 

207 

0 

.792 

2.844 

209 

0 

.813 

2.969 

126 

0 

.938 

2. 9B666666 

178 

0 

.938 

2.96 

208 

0 

.813 

2.844 

206 

0 

.771 

2.844 

50 

0 

1.0 

3.0 

6 

0 

1.0 

2.782 

75 

0 

1.0 

1.0 

amATERIAL  PROPERTIES  ARE  DEFINED 
MATERIAL  1  THERMO -ELASTIC 
-2000  29600000  .3  .0000065 

2000  29600000  .3  .0000065 

ATHE  element  is  defined 

EGROUP  1  PLANE  5U-AXI5YMM 

ATHE  MESH  FOR  THE  TUBE  IS  GENERATED 

GSURFACE  3  1  37  39  ELI-1  EL2-3  NO-9  NC-B 

GSURFACE  39  37  60  61  ELI-1  EL2-2  NO-9  NC=B 

GSURFACE  61  60  70  71  ELI-1  EL2-4  NO-9  NC=B 

GSURFACE  71  70  B5  B7  ELI-1  EL2-1  NO-9  NC-B 

ANODES  FOR  THE  WELD  ARE  GENERATED 

LINE  A  3  100  NCENTER-200  EL-?  M-l 

LINE  N  100  204  126  152  178 

LINE  A  204  211  N CENTER- 209  EL-1  M-l 

LINE  C  100  211  204 

LINE  A  211  210  NCENTER-209  EL-1  M-l 
LINE  N  210  39  20B  207  206 
LINE  C  211  39  210 

ATHE  MESH  FOR  THE  WELD  IS  GENERATED 

GSURFACE  100  3  39  211  ELI-3  EL2-3  NO-9  NC-B 

*NODES  FOR  THE  HEADER  ARE  GENERATED 

LINE  C  100  39  204  211  210 

LINE  S  50  6  EL-3  M-l 

LINE  S  6  61  EL-3  M-l 

LINE  C  50  61  6 

*THE  MESH  FOR  THE  HEADER  IS  GENERATED 
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GSURFACE  50  100  39  61  EL1*2  EL2*6  NO»9  NC=B 
GSURFACE  6  61  72  75  EL1  =  3  EL2*4  N0*9  NC=N 
^PORTIONS  OF  THE  MODEL  ARE  DEFINED  FOR  PLOTTING 
EZONE  TUBE01 

1  STEP  1  TO  3 
EZONE  TUBE02 

4  STEP  1  TO  10 
EZONE  HELD 

11  STEP  1  TO  19 
EZONE  HEAD 

20  STEP  1  TO  43 

BZONE  CLOSE  XMIN=0  XMAX=0  YMIN=.63  YMAX=1.0  ZMIN=2.78~  ZMAX=3.1 
*THE  MASTER  DEGREES  OF  FREEDOM  ARE  DEFINED 
MASTER  REACTI ONS  = YES  IDOF* 100111  DT=1  NSTEP=2  ITP56  =  1 
PRINTOUT  MAX 
ANALYSIS  TYPE=STATIC 
aTHE  LOAD  IS  APPLIED  IN  2 'STEPS 
TIMEFUNCTION 
0  0 

2  1 
A 

TOL  T=E  1*15  DN*. 00001 

*NODES  ONLY  USED  FOR  GENERATION  ARE  FIXED,  THE  BOTTOM  RIGHT  HAND 

*NODE  OF  THE  HEADER  IS  FIXED 

BOUNDARIES  111111  TYPE=NODES  /200  209  75 

a THE  NODES  ON  THE  RIGHT  HAND  SIDE  OF  THE  HEADER  ARE  CONSTRAINED  TO 
AROLL  ONLY  IN  THE  Z  DIRECTION 

BOUNDARIES  110111  TYPE = NODES  /  50  101  STEP  1  TO  105  6 
BOUNDARIES  110111  TYPE’NODES  /  148  156  163  170  177  1B5  192 
*THE  CONTACT  SURFACE  IS  DEFINED 
C GROUP  1  C0NTACT2  S=AXISYM 
CONTACTSURFACE  1 

61  155  162  169  176  184  191  198  72 
CONTACTSURFACE  2 

87  56  71  52  48  45  42  38  34  31  61 
CONTACTPAIR  1  T=2  C=>1  F=0 
aTHE  INTERNAL  PRESSURE  LOADING 
LOADS  ELEMENT 
1  -1  700 
STEP  1  TO 

10  -1  700 

1  2  700 

13  2  700 

12  2  700 

11  2  700 

21  2  700 
20  2  700 

* THE  LONGITUDINAL  TUBE  LOAD 

LOADS  CONCENTRATED 

85  3  -25.3 

59  3  -101.3 

87  3  -25.3 

ATHE  THERMAL  LOAD 

LOADS  TEMP  TR=0 

50  350.0 

101  348.1833333 

102  346.3666667 

103  344.55 

104  342.7333333 

105  340.9166667 


Si 


sk«£ as 


6 

339.1 

148 

327.9625 

156 

316.325 

163 

305.6875 

170 

294.55 

177 

283.4125 

185 

272.275 

192 

261.1375 

75 

1 

350 

7 

350 

10 

350 

13 

350 

16 

350 

19 

350 

37 

350 

23 

350 

26 

350 

29 

350 

60 

350 

33 

350 

36 

350 

41 

350 

44 

350 

47 

350 

51 

350 

54 

350 

70 

350 

58 

350 

85 

350 

2 

350 

3 

350 

62 

350 

63 

350 

64 

350 

65 

350 

66 

350 

100 

350 

113 

350 

112 

350 

111 

350 

PA  TAR  CREA3H 

* INPUT  FOR  THE  FINE  MESH,  2-DIMENSIONAL  AXISYMMETRIC  MODEL 
HEAD  'HEADER' 

CONTROL  B-4800  SYSTEM  => 2  DEVICE»2  OPTION-1024 
CONTROL  PLOTE=WAIT  PLOTS=NO 
^COORDINATES  ARE  DEFINED 

COORDINATES  /  ENTRIES  NODE  X  Y  Z 


1 

0 

.63 

3.0 

3 

0 

.693 

3.0 

37 

0 

.63 

2.844 

39 

0 

.75 

2.844 

60 

0 

.63 

2.782 

61 

0 

.75 

2.782 

70 

0 

.63 

1.0 

71 

0 

.75 

1.0 

72 

0 

.75001 

1.0 

85 

0 

.63 

0.0 

87 

0 

.75 

0.0 

100 

0 

.938 

3.0 

200 

0 

.82 

2.837 

152 

0 

.938 

2.97333333 

256 

0 

.901 

2.88 

248 

0 

.792 

2.844 

78 


251 

0 

.813 

2.969 

126 

0 

.938 

2.98666666 

178 

0 

.938 

2.96 

250 

0 

.813 

2.844 

246 

0 

.771 

2.844 

50 

0 

1.0 

3.0 

6 

0 

1.0 

2.782 

75 

0 

1.0 

1.0 

113 

0 

.938 

2.99333333 

139 

0 

.938 

2.98 

165 

0 

.938 

2.96666666 

249 

0 

.8025 

2.844 

247 

0 

.7815 

2.844 

245 

0 

.7605 

2.844 

★MATERIAL  PROPERTIES  ARE  DEFINED 
MATERIAL  1  ELASTIC  E-29600000  NU-.3 
★THE  ELEMENT  IS  DEFINED 
EGROUP  1  PLANE  SU-AXISYMM 
★THE  MESH  FOR  THE  TUBE  IS  GENERATED 
GSURFACE  3  1  37  39  EL1=2  EL2-6  NO-9 
GSURFACE  39  37  60  61  ELI-2  EL2-4  NO-9 

GSURFACE  61  60  70  71  EL1=2  EL2-8  NO-9 

GSURFACE  71  70  35  37  ELI-2  EL2-2  NO-9 

★NODES  FOR  THE  WELD  ARE  GENERATED 
LINE  A  3  100  NCENTER-200  EL-6  M-l 
LINE  N  100  178  113  126  139  15.2  165 
LINE  A  178  256  NCENTER-251  EL=3  M-l 
DINE  C  100  256  178 
LINE  A  256  250  NCENTER-251  EL-3  M-l 
LINE  N  250  39  249  248  247  246  245 
LINE  C  256  39  250 

★THE  MESH  FOR  THE  WELD  IS  GENERATED 
GSURFACE  100  3  39  256  ELI-6  EL2=6  NO=9 
★NODES  FOR  THE  HEADER  ARE  GENERATED 
LINE  C  100  39  178  256  250 
LINE  S  50  6  EL=6  M=1 
LINE  S  6  61  EL=6  M-l 

LINE  C  SO  61  6 

★THE  MESH  FOR  THE  HEADER  IS  GENERATED 
GSURFACE  50  100  39  61  ELI-4  EL2-12  NO-9 
GSURFACE  6  61  72  75  ELI -6  EL2-8  NO-9  NC-N 
★PORTIONS  OF  THE  TUBE  ARE  DEFINED  FOR-PtiOTgHNG 
EZONE  TUBE01 

1  STEP  1  TO  12 
EZONE  TUBE02 

13  STEP  1  TO  44 
EZONE  WELD 
45  STEP  1  TO  80 
EZONE  HEAD 
81  yPE?P  l  TO  256 

BZONE  CLOSE  XMIN-0  XMAX-0  YMIN-.63  YMAX- 1.094  ZMIN-2.688  ZMAX-3.1 
BZONE  HEADOl  XMIN-0  XMAX-0  YMIN-1.094  YMAX-2. 0  ZMIN-2.688  ZMAX-3.0 

BZONE  HEAD02  xnIN-0  XMAX-0  YMIN-1.094  YMAX-2. 0  ZMIN-1.0  ZMAX-2.688 

BZONE  HEAD03  XMIN-0  XMAX-0  YMIN-.63  YMAX-1.094  ZMIN-1.0  ZMAX-2.688 

BZONE  TUBE  XMIN-0  XMAX-0  YMIN- . 63  YMAX-.693  ZMIN-0  ZMAX-1.0 
★THE  MASTER  DEGREES  OF  FREEDOM  ARE  DEFINED 
MASTER  REACTIONS -YES  ID0F-100111  DT-1  NSTEP-2 
PRINTOUT  MAX 
ANALYSIS  TYPE-STATIC 
★THE  LOAD  IS  APPLIED  IN  2  STEPS 
TIMEFUNCTION 
0  0 

2  1 
•M 

TOL  T-E  1=15  DN-. 00001 

★NODES  USED  FOR  GENERATION  ARE  FIXED,  THE  BOTTOM  RIGHT  HAND  NODE 
★OF  THE  HEADER  IS  FIXED 

BOUNDARIES  111111  TYPE-NODES  /  200  251  75 
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*THE  NODES  AT  THE  RIGHT  HAND  SIDE  OF  THE  HEADER  ARE  CONSTRAINED  TO 
*ROLL  ONLY  IN  THE  Z -DIRECTI ON 

BOUNDARIES  110111  TYPE-NODES/  50  368  STEP  1  TO  378  6  558  STEP  13  TO  740 
*THE  CONTACT  SURFACES  ARE  DEFINED 
CGROUP  1  C0NTACT2  S-AXISYM 
CONTACTSURFACE  1 

87  210  205  199  71  191  186  181  175  170  164  159  154  148  143  137  132  127  121  116  61 
CONTACTSURFACE  2 

61  570  583  596  609  622  635  648  661  674  687  700  713  726  739  752  72 

CONTACTPAIR  1  T*2  C«1  F«0 

*THE  INTERNAL  PRESSURE  LOAD 

LOADS  ELEMENT 

2  -1  700 

STEP  2  TO 

40  -1  700 

2  2  700 

1  2  700 

41  2  700 
STEP  1  TO 
46  2  700 
77  2  700 
STEP  1  TO 
80  2  700 

*THE  LONGITUDINAL  TUBE  LOAD 
LOADS  CONCENTRATED 
85  3  -12.66666666 

217  3  -50.66666666 

216  3  -25.33333333 

215  3  -50.66666666 

87  3  -12.66666666 

*THE  THERMAL  LOAD 
LOADS  TEMP  TR=0 
50  350 

368  349.0916667 

369  348.1833333 

370  347.275 

371  346.3666667 

372  345.4583333 

373  344.55 

374  343.6416667 

375  342.7333333 

376  341.825 

377  340.9166667 

378  340.0083333 

6  339.1 

558  333.53125 

571  327.9625 

584  322.39375 

597  316.825 

610  311.25625 

623  305.6875 

636  300.11875 

649  294.55 

662  288.98125 

675  283.4125 

688  277.84375 

701  272.275 

714  266.70625 

727  261.1375 

740  255.56875 

75  250 

1  350 

STEP  1  TO 


STEP  1  TO 


APPENDIX  H.  ADINA-IN  INPUT  FILE  FOR  THE  3-DIMENSIONAL 

CONTACT  SURFACE 


’"'THESE  ARE  THE  CONTACT  SURFACES  FOR  THE  3 -DIMENSIONAL  MODEL 
CGROUP  1  CONTACT 3 

’"'CONTACTSURFACE  1  IS  THE  HEADER  AROUND  THE  INNER  TUBE 
CONTACTSURFACE  1 


1 

123 

124 

970 

969 

2 

969 

970 

982 

981 

3 

981 

982 

994 

993 

4 

993 

994 

69 

68 

5 

124 

125 

971 

970 

6 

970 

971 

983 

982 

7 

982 

983 

995 

994 

8 

994 

995 

70 

69 

9 

125 

126 

1021 

971 

10 

971 

1021 

1029 

983 

11 

983 

1029 

1037 

995 

12 

995 

1037 

71 

70 

13 

126 

127 

1022 

1021 

14 

1021 

1022 

1030 

1029 

15 

1029 

1030 

1038 

1037 

16 

1037 

1038 

72 

71 

17 

127 

128 

1057 

1022 

18 

1022 

1057 

1065 

1030 

19 

1030 

1065 

1073 

1038 

20 

1038 

1073 

73 

72 

21 

128 

129 

1058 

1057 

22 

1057 

1058 

1066 

1065 

23 

1065 

1066 

1074 

1073 

24 

1073 

1074 

74 

73 

25 

129 

130 

1093 

1058 

26 

1058 

1093 

1101 

1066 

27 

1066 

1101 

1109 

1074 

28 

1074 

1109 

75 

74 

29 

130 

131 

1094 

1093 

30 

1093 

1094 

1102 

1101 

31 

1101 

1102 

1110 

1109 

32 

1109 

1110 

76 

75 

33 

131 

132 

1129 

1094 

34 

1094 

1129 

1137 

1102 

35 

1102 

1137 

1145 

1110 

36 

1110 

1145 

77 

76 

37 

132 

11 

1130 

1129 

38 

1129 

1130 

1138 

1137 

39 

1137 

1138 

1146 

1145 

40 

1145 

1146 

6 

77 

82 


■ 


41  11  122  1161 

42  1130  1161  1165 

43  1138  1165  1169 

44  1146  1169  67 

45  122  123  969 

46  1161  969  981 

47  1165  981  993 

48  1169  993  68 

*CONTACTSURFACE  2 
CONTACTSURFACE  2 

1  585  586  1354 

2  1323  1354  1358 

3  1331  1358  1362 

4  1339  1362  531 

5  586  587  1180 

6  1354  1180  1188 

7  1358  1188  1196 

8  1362  1196  532 

9  587  588  1181 

10  1180  1181  1189 

11  1188  1189  1197 

12  1196  1197  533 

13  588  589  1182 

14  1181  1182  1190 

15  1189  1190  1198 

16  1197  1198  534 

17  589  590  1212 

18  1182  1212  1218 

19  1190  1218  1224 

20  1198  1224  535 

21  590  591  1213 

22  1212  1213  1219 

23  1218  1219  1225 

24  1224  1225  536 

25  591  592  1250 

26  1213  1250  1258 

27  1219  1258  1266 

28  1225  1266  537 

29  592  593  1251 

30  1250  1251  1259 

31  1258  1259  1267 

32  1266  1267  538 

33  593  594  1286 

34  1251  1286  1294 

35  1259  1294  1302 

36  1267  1302  539 

37  594  473  1287 

38  1286  1287  1295 

39  1294  1295  1303  ; 

40  1302  1303  468 

41  473  584  1322  ! 


L  1130 
5  1138 
)  1146 
’  6 
>  1161 
.  1165 
I  1169 
t  67 

:  IS  THE  HEADER  AROUND  THE  OUTER  TUBE 

1323 

1331 

1339 

530 
1354 
1358 
1362 

531 
1180 
1188 

1196 

532 
1181 

1189 

1197 

533 
1182 

1190 

1198 

534 
1212 
1218 

1224 

535 
1213 
1219 

1225 

536 

1250 

1258 
1266 

537 

1251 

1259 
1267 

538 
1286 
1294 
1302 

539 
1287 


42  1287  1322  1330  1295 

43  1295  1330  1338  1303 

44  1303  1338  529  468 

45  584  585  1323  1322 

46  1322  1323  1331  1330 

47  1330  1331  1339  1338 

48  1338  1339  530  529 

*CONTACTSURFACE  3  IS  THE  INNER  TUBE 
CONTACTSURFACE  3 

1  195  237  238  196 

2  237  90  91  238 

3  90  297  298  91 

4  297  339  340  298 

5  339  381  382  340 

6  381  68  69  382 

7  196  238  239  197 

8  238  91  92  239 

9  91  298  299  92 

10  298  340  341  299 

11  340  382  383  341 

12  382  69  70  383 

13  197  239  240  198 

14  239  92  93  240 

15  92  299  300  93 

16  299  341  342  300 

17  341  383  384  342 

18  383  70  71  384 

19  198  240  241  199 

20  240  93  94  241 

21  93  300  301  94 

22  300  342  343  301 

23  342  384  385  343 

24  384  71  72  385 

25  199  241  242  200 

26  241  94  95  242 

27  94  301  302  95 

28  301  343  344  302 

29  343  385  386  344 

30  385  72  73  386 

31  200  242  243  201 

32  242  95  96  243 

33  95  302  303  96 

34  302  344  345  303 

35  344  386  387  345 

36  386  73  74  387 

37  201  243  244  202 

38  243  96  97  244 

39  96  303  304  97 

40  303  345  346  304 

41  345  387  388  346 

42  387  74  75  388 

43  202  244  245  203 


44  244  97  98  245 

45  97  304  305  98 

46  304  346  347  305 

47  346  388  389  347 

48  388  75  76  389 

49  203  245  246  204 

50  245  98  99  246 

51  98  305  306  99 

52  305  347  348  306 

53  347  389  390  348 

54  389  76  77  390 

55  204  246  235  193 

56  246  99  8  235 

57  99  306  295  8 

58  306  348  337  295 

59  348  390  379  337 

60  390  77  6  379 

61  193  235  236  194 

62  235  8  89  236 

63  8  295  296  89 

64  295  337  338  296 

65  337  379  380  338 

66  379  6  67  380 

67  194  236  237  195 

68  236  89  90  237 

69  89  296  297  90 

70  296  338  339  297 

71  338  380  381  339 

72  380  67  68  381 

*  CONTACT  SURF  ACE  4  IS  THE  OUTER  TUBE 
CONTACTSURFACE  4 

1  657  699  700  658 

2  699  552  553  700 

3  552  759  760  553 

4  759  801  802  760 

5  801  843  844  802 

6  843  530  531  844 

7  658  700  701  659 

8  700  553  554  701 

9  553  760  761  554 

10  760  802  803  761 

11  802  844  845  803 

12  844  531  532  845 

13  659  701  702  660 

14  701  554  555  702 

15  554  761  762  555 

16  761  803  804  762 

17  803  845  846  804 

18  845  532  533  846 

19  660  702  703  661 

20  702  555  556  703 

21  555  762  763  556 


22  762  804  805  763 

23  804  846  847  805 

24  846  533  534  847 

25  661  703  704  662 

26  703  556  557  704 

27  556  763  764  557 

28  763  805  806  764 

29  805  847  848  806 

30  847  534  535  848 

31  662  704  705  663 

32  704  557  558  705 

33  557  764  765  558 

34  764  806  807  765 

35  806  848  849  807 

36  848  535  536  849 

37  663  705  706  664 

38  705  558  559  706 

39  558  765  766  559 

40  765  807  808  766 

41  807  849  850  808 

42  849  536  537  850 

43  664  706  707  665 

44  706  559  560  707 

45  559  766  767  560 

46  766  808  809  767 

47  808  850  851  809 

48  850  537  538  .851 

49  665  707  708  566 

50  707  560  561  708 

51  560  767  768  561 

52  767  809  810  768 

53  809  851  852  810 

54  851  538  539  852 

55  666  708  697  655 

56  708  561  470  697 

57  561  768  757  470 

58  76 Q  81 0  799  757 

59  810  852  841  799 

60  852  539  468  841 

61  655  697  698  656 

62  697  470  551  698 

63  470  757  758  551 

64  757  799  800  758 

65  799  841  842  800 

66  841  468  529  842 

67  656  698  699  657 

68  698  551  552  699 

69  551  758  759  552 

70  758  800  801  759 

71  800  842  843  801 

72  842  529  530  843 
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APPENDIX  I.  ADINA-IN  INPUT  FILE  FOR  THE  3-DIMENSIONAL 

PROBLEM 


DATAB  CREATE 
HEAD  '3D-HEADER' 

CONTROL  B=4800  SYSTEM=2  DEVICE=2  0PTION*1024 
CONTROL  PLOTE=HAXT  PLOTS=NO 

*THE  COORDINATES  FOR  THE  INNER  TUBE  ARE  DEFINED 
SYSTEM  1  TYF£=CYLINDRICAL  X=9.4  Y=1  2=1.125 


COORDINATES  /ENTRIES 

NODE 

R 

THETA 

XL 

1 

.63 

0 

0 

2 

.75 

0 

0 

3 

.63 

0 

.156 

4 

.75 

0 

.156 

5 

.63 

0 

.218 

6 

.75 

0 

.218 

7 

.63 

0 

2.0 

8 

.75 

0 

2.0 

9 

.63 

0 

4.0 

10 

.75 

0 

4.0 

11 

.75001 

0 

2.0 

LINE 

CYLINDRICAL 

1 

1 

6  1 

LINE 

CYLINDRICAL 

2 

2 

6  1 

LINE 

CYLINDRICAL 

3 

3 

6  1 

LINE 

CYLINDRICAL 

4 

4 

6  1 

LINE 

CYLINDRICAL 

5 

5 

6  1 

LINE 

CYLINDRICAL 

6 

6 

6  1 

. 

LINE 

CYLINDRICAL 

7 

7 

6  1 

LINE 

CYLINDRICAL 

8 

a 

6  1 

LINE 

CYLINDRICAL 

9 

9 

6  1 

LINE 

CYLINDRICAL 

10 

10  6  1 

LINE 

CYLINDRICAL 

11 

11  6  1 

*THE 

MASTER  DEGREES 

OF  FREEDOM  ARE 

DEFINED 

MASTER  REACTIONS = YES  IDOF-OOOlll 

PRINTOUT  V=MIN  10=0 

ANALYSIS  TYPE = STATIC 

VIEW  1  X=1 

VIEW  2  X=-l 

VIEW  3  Y=1 

VIEW  4  Y=-l 

VIEW  5  2=1 

VIEW  6  2  —  1 

*THE  MATERIAL  PROPERTIES  ARE  DEFINED 
MATERIAL  1  ELASTIC  E=29600000  NU=.3 
*THE  ELEMENT  IS  DEFINED 
EGROUP  1  SOLID 

STRESSTABLE  112345678  21 
EDATA/ ENTRIES  EL  TABLE  PRINT 

1  1  YES 

STEP  1  TO 
288  1  YES 

*THE  MESH  FOR  THE  INNER  TUBE  IS  GENERATED 
GVOLUME  10  10  9  9  8  8  7  7  EL1=6  EL2=1  EL3=4  N0=20 
GVOLUME  88776655  EL1=6  EL2=1  EL3=4  N0=20 
GVOLUME  66554433  EL1=6  EL2=1  EL3=1  N0=20 
GVOLUME  44332211  EL1=6  EL2=1  EL3=1  N0=20 
*THE  COORDINATES  FOR  THE  OUTER  TUBE  ARE  DEFINED 
SYSTEM  2  TYPE =CYL INDR I CAL  X=9.4  Y=3.0  Z  =  1.125 
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COORDINATES  /ENTRIES  NODE  R  THETA  XL 

463  .63  0  -.4564 

464  .75  0  -.49 

465  .63  0  -.09327705 

466  .75  0  -.11162873 


479 

.63 

180 

-.1036 

490 

.75 

180 

-.07 

467 

.63 

0 

.218 

468 

.75 

0 

.218 

469 

.63 

0 

1.64381482 

470 

.75 

0 

1.61759259 

471 

.63 

0 

4.0 

472 

.75 

0 

4.0 

473 

.75001 

0 

1.61759041 

589 

.75001 

180 

1.94537256 

474 

.63 

30 

-.43276688 

475 

.63 

60 

-.3682 

476 

.63 

90 

-.28 

477 

.63 

120 

-.1918 

478 

.63 

150 

-.12723312 

480 

.63 

210 

-.12723312 

481 

.63 

240 

-.1918 

482 

.63 

270 

-.28 

483 

.63 

300 

-.3682 

484 

.63 

330 

-.43276688 

485 

.75 

30 

-.46186535 

486 

.75 

60 

-.385 

487 

.75 

90 

-.28 

488 

.75 

120 

-.175 

489 

.75 

150 

-.09813465 

491 

.75 

210 

-.09813465 

492 

.75 

240 

-.175 

493 

.75 

270 

-.28 

494 

.75 

300 

-.385 

495 

.75 

330 

-.46186535 

584 

.75001 

30 

1.63954765 

585 

.75001 

60 

1.69953594 

586 

.75001 

90 

1.78148148 

587 

.75001 

120 

1.86342702 

588 

.75001 

150 

1.92341531 

590 

.75001 

210 

1.92341531 

591 

.75001 

240 

1.86342702 

592 

.75001 

270 

1.78148148 

593 

.75001 

300 

1.69953594 

594 

.75001 

330 

1.63954765 

496 

.63 

30 

-.08036908 

497 

.63 

60 

-.04510387 

498 

.63 

90 

.0030693 

499 

.63 

120 

.05124247 

500 

.63 

150 

.08650768 

501 

.63 

180 

.09941564 

502 

.63 

210 

.08650768 

503 

.63 

240 

.05124247 

504 

.63 

270 

.0030693 

505 

.63 

300 

-.04510387 

506 

.63 

330 

-.08036908 

507 

.75 

30 

-.09626212 

508 

.75 

60 

-.05427972 

509 

.75 

90 

.0030693 

510 

.75 

120 

.06041831 

511 

.75 

150 

.10240071 

512 

.75 

180 

.11776732 

513 

.75 

210 

.10240071 

514 

.75 

240 

.06041831 

515 

.75 

270 

. 0030693 

516 

.75 

300 

-.05427972 

517 

.75 

330 

-.09626212 

540 

.63 

30 

1.66225867 

541 

.63 

60 

1.71264813 

542 

.63 

90 

1.78148148 

543 

.63 

120 

1.85031484 

544 

.63 

150 

1.9007043 

545 

.63 

180 

1.91914817 

546 

.63 

210 

1.9007043 

547 

.63 

240 

1.85031484 

548 

.63 

270 

1.78148148 

549 

.63 

300 

1.71264813 

550 

.63 

330 

1.66225867 

551 

.75 

30 

1.63954951 

552 

.75 

60 

1.69953704 

553 

.75 

90 

1.78148148 

554 

.75 

120 

1.86342593 

555 

.75 

150 

1.92341346 

556 

.75 

180 

1.94537037 

557 

.75 

210 

1.92341346 

558 

.75 

240 

1.86342593 

559 

.75 

270 

1.78148148 

560 

.75 

300 

1.69953704 

561 

.75 

330 

1.63954951 

LINE  NODES  463  463  474  475 

476  STEP  1  TO  484 
LINE  NODES  464  464  485  486 

487  STEP  1  TO  495 
LINE  NODES  465  465  496  497 

498  STEP  1  TO  506 
LINE  NODES  466  466  507  508 

509  STEP  1  TO  517 
LINE  CYLINDRICAL  467  467  6  1 
LINE  CYLINDRICAL  468  468  6  1 
LINE  NODES  469  469  540  541 

542  STEP  1  TO  550 
LINE  NODES  470  470  551  552 

553  STEF  1  TO  561 
LINE  CYLINDRICAL  471  471  6  1 
LINE  CYLINDRICAL  472  472  6  1 
LINE  NODES  473  473  584  585 

586  STEP  1  TO  594 

*THE  MESH  FOR  THE  OUTER  TUBE  IS  GENERATED 

GVOLUME  472  472  471  471  470  470  469  469  ELI-6  EL2-1  EL3-4  NO-20 

GVOLUME  470  470  469  469  468  468  467  467  EL1*6  EL2-1  EL3-4  N0*20 

GVOLUME  468  468  467  467  466  466  465  465  EL1*6  EL2“1  EL3  =  1  N0*20 

GVOLUME  466  466  465  465  464  464  463  463  EL1*6  EL2-1  EL3*1  N0*20 

*THE  COORDINATES  FOR  THE  HEADER  AROUND  THE  TUBES  ARE  DEFINED 
SYSTEM  N=3  TYPE*CARTE2 
COORDINATES  /  ENTRIES 


AN  X* 

9.4  Y  =  0 

2*0 

NODE 

5 

X 

Y 

z 

925 

0.0 

0.0 

0.0 

926 

.156 

0.0 

0.0 

927 

.218 

0.0 

0.0 

928 

2.0 

0.0 

0.0 

929 

0.0 

0.0 

1.125 

930 

.156 

0.0 

1.125 

931 

.218 

0.0 

1.125 

932 

2.0 

0.0 

1.125 

933 

0.0 

0.0 

2.25 

934 

.156 

0.0 

2.25 

935 

.218 

0.0 

2.25 

936 

2.0 

0.0 

2.25 

937 

0.0 

2.0 

0.0 

938 

.156 

2.0 

0.0 

939 

.218 

2.0 

0.0 

940 

2.0 

2.0 

0.0 

941 

0.0 

2.0 

1.125 

942 

.156 

2.0 

1.125 

943 

.218 

2.0 

1.125 

944 

2.0 

2.0 

1.125 

945 

0.0 

2.0 

2.25 

946 

.156 

2.0 

2.25 

947 

.218 

2.0 

2.25 

948 

2.0 

2.0 

2.25 

*NOD£S  OF  THE  CENTERS  OF  THE  TUBES  FOR  MESH  GENERATION 


* OUTER  TUBE  HEADER  NODES 


949 

0.0 

1.0 

1.125 

950 

.156 

1.0 

1.125 

951 

.218 

1.0 

1.125 

952 

2.0 

1.0 

1.125 

953 

0.0 

3.0 

1.125 

954 

.156 

3.0 

1.125 

955 

.218 

3.0 

1.125 

956 

NODES 

2.0 

3.0 

1.125 

1230 

1.41 

4.7 

2.25 

1231 

.218 

4.2764467 

2.25 

1232 

-.171 

4.13822335 

2.25 

1233 

-.56 

4.0 

2.25 

1234 

1.41 

4.7 

1.125 

1235 

.218 

4.2764467 

1.125 

1236 

-.171 

4.13822335 

1.125 

1237 

-.56 

4.0 

1.125 

1238 

1.41 

4.7 

0.0 

1239 

.218 

4.2764467 

0.0 

1240 

-.171 

4.13822335 

0.0 

1241 

-.56 

4.0 

0.0 

MESH  FOR 
124 

THE  INNER  TUBE  HEADER 

LINE  NODES  123  125  124 
LINE  ARC  68  70  NCEN-951  EL-1  M-l  NCO-ALL 

LINE  ARC  46  48  NCEN-950  EL-1  M-l  NCO-ALL 

LINE  ARC  24  26  NCEN-949  EL-1  M-l  NCO-ALL 

GVOLUME  948  936  125  123  947  935  70  68  ELI-1  EL2-1  EL3-4  NC-N 

GVOLUME  947  935  70  68  946  934  48  46  ELI-1  EL2-1  EL3-1  NC-A 

GVOLUME  946  934  48  46  945  933  26  24  0.1-1  EL2-1  EL3-1  NC-A 

* 

LINE  NODES  125  127  126 

LINE  ARC  70  72  NCEN-951  EL-1  M-l  NCO-ALL 

LINE  ARC  48  50  NCEN-950  EL-1  M-l  NCO-ALL 

LINE  ARC  26  28  NCEN-949  EL-1  M-l  NCO-ALL 

GVOLUME  936  932  127  125  935  931  72  70  ELI-1  EL2-1  EL3-4  NC-N 

GVOLUME  935  931  72  70  934  930  50  48  ELI-1  EL2-1  EL3-1  NC-A 

GVOLUME  934  930  50  48  933  929  28  26  ELI-1  EL2-1  EL3-1  NC-A 


LINE  NODES  127  129  128 

LINE  ARC  72  74  NCEN-951  EL-1  M-l  NCO-ALL 

LINE  ARC  50  52  NCEN-950  EL-1  M-l  NCO-ALL 

LINE  ARC  28  30  NCEN-949  EL-1  M-l  NCO-ALL 

GVOLUME  932  928  129  127  931  927  74  72  ELI-1  EL2-1  EL3-4  NC-N 
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GVOLUME  931  927  74  72  930  926  52  50  ELI-1  EL2-1  EL3-1  NC-A 

GVOLUME  930  926  52  50  929  925  30  28  ELI-1  EL2-1  EL3-1  NC-A 

a 

LINE  NODES  129  131  130 

LINE  ARC  74  76  NCEN-951  EL-1  M-l  NCO-ALL 

LINE  ARC  52  54  NCEN-950  EL-1  M-l  NCO-ALL 

LINE  ARC  30  32  NCEN-949  EL-1  M-l  NCO-ALL 

GVOLUME  928  940  131  129  927  939  76  74  ELI-1  EL2-1  EL3-4  NC-N 

GVOLUME  927  939  76  74  926  938  54  52  ELI-1  EL2-1  EL3-1  NC-A 

GVOLUME  926  938  54  52  925  937  32  30  ELI-1  EL2-1  EL3-1  NC-A 

A 

LINE  NODES  131  11  132 

LINE  ARC  76  6  NCEN-951  EL-1  M-l  NCO-ALL 

LINE  ARC  54  4  NCEN-950  EL-1  M-l  NCO-ALL 

LINE  ARC  32  2  NCEN-949  EL-1  M-l  NCO-ALL 

GVOLUME  940  944  11  131  939  943  6  76  ELI-1  EL2-1  EL3-4  NC-N 

GVOLUME  939  943  6  76  938  942  4  54  ELI-1  EL2-1  EL3-1  NC-A 

GVOLUME  938  942  4  54  937  941  2  32  ELI-1  EL2-1  EL3-1  NC-A 

A 

LINE  NODES  11  123  122 

LINE  ARC  6  68  NCEN-951  EL-1  M-l  NCO-ALL 

LINE  ARC  4  46  NCEN-950  EL-1  M-l  NCO-ALL 

LINE  ARC  2  24  NCEN-949  EL-1  M-l  NCO-ALL 

GVOLUME  944  948  123  11  943  947  68  6  ELI-1  EL2-1  EL3-4  NC-N 
GVOLUME  943  947  68  6  942  946  46  4  ELI-1  EL2-1  EL3-1  NC-A 

GVOLUME  942  946  46  4  941  945  24  2  ELI-1  EL2-1  EL3-1  NC-A 

* GENERATION  OF  THE  MESH  FOR  THE  OUTER  TUBE  HEADER 
LINE  NODES  587  589  588 
LINE  NGDES  532  534  533 
LINE  NODES  510  512  511 
LINE  NODES  488  490  489 

GVOLUME  948  944  589  587  947  943  534  532  ELI-1  EL2-1  EL3-4  NC-N 
LINE  NODES  532  510  891 

GVOLUME  947  943  534  532  946  942  512  510  ELI-1  EL2-1  EL3-1  NC-A 
LINE  NODES  510  488  909 

GVOLUME  946  942  512  510  945  941  490  488  ELI-1  EL2-1  EL3-1  NC-A 

A 

LINE  NODES  589  591  590 
LINE  NODES  534  536  535 
LINE  NODES  512  514  513 
LINE  NODES  490  492  491 

GVOLUME  944  940  591  589  943  939  536  534  ELI-1  EL2-1  EL3-4  NC-N 
LINE  NODES  536  514  893 

GVOLUME  943  939  536  534  942  938  514  512  ELI-1  EL2-1  EL3-1  NC-A 
LINE  NODES  514  492  911 

GVOLUME  942  938  514  512  941  937  492  490  ELI-1  EL2-1  EL3-1  NC-A 

A 

LINE  NODES  591  593  592 
LINE  NODES  536  538  537 
LINE  NODES  514  516  515 
LINE  NODES  492  494  493 

GVOLUME  940  1238  593  591  939  1239  538  536  ELI-1  EL2-1  EL3-4  NC-N 
LINE  NODES  538  516  894 

GVOLUME  939  1239  538  536  938  1240  516  514  ELI-1  EL2-1  EL3-1  NC-A 
LINE  NODES  516  494  912 

GVOLUME  938  1240  516  514  937  1241  494  492  ELI-1  EL2-1  EL3-1  NC-A 

A 

LINE  NODES  593  473  594 
LINE  NODES -538  468  539 
LINE  NODES'  516  466  517 


LINE  NODES  494  464  495 

GVOLUME  1238  1234  473  593  1239  1235  468  538  ELI-1  EL2»1  EL3-4  NON 

GVOLUME  1239  1235  468  538  1240  1236  466  516  EL1=1  EL2-1  EL3-1  NOA 

GVOLUME  1240  1236  466  516  1241  1237  464  494  EL1=1  EL2»1  EL3«1  NOA 

* 

LINE  NODES  473  585  584 
LINE  NODES  468  530  529 
LINE  NODES  466  508  507 
LINE  NODES  464  486  485 

GVOLUME  1234  1230  585  473  1235  1231  530  468  ELI-1  EL2-1  EL3*4  NON 
LINE  NODES  530  508  890 

GVOLUME  1235  1231  530  468  1236  1232  508  466  EL1*1  EL2-1  EL3=1  NC=>A 
LINE  NODES  508  486  908 

GVOLUME  1236  1232  508  466  1237  1233  486  464  EL1  =  1  EL2-1  EL3»1  NOA 
* 

LINE  NODES  585  587  586 
LINE  NODES  530  532  531 
LINE  NODES  508  510  509 
LINE  NODES  486  488  487 

GVOLUME  1230  948  587  585  1231  947  532  530  EL1»1  EL2*1  EL3-4  NON 

GVOLUME  1231  947  532  530  1232  946  510  508  ELI-1  EL2*1  EL3=*1  NOA 

GVOLUME  1232  946  510  508  1233  945  488  486  ELI-1  EL2-1  EL3*1  NOA 

★COORDINATES  FOR  THE  REST  OF  THE  HEADER  ARE  DEFINED 

SYSTEM  N=0 


COORDINATES/  ENTRIES  NODE 

X 

Y 

Z 

1366 

9.09 

6.32 

0 

1367 

9.09 

6.32 

2.25 

1368 

7.86 

4.67 

0 

1369 

7.86 

4.67 

2.25 

1370 

6.9 

6.79 

0 

1371 

6.9 

6.79 

2.25 

1372 

6.644 

4.85 

0 

1373 

6.644 

4.85 

2.25 

1374 

4.91 

6.37 

0 

1375 

4.91 

6.37 

2.25 

1376 

5.39 

4.6 

0 

1377 

5.39 

4.6 

2.25 

1378 

3.16 

5.575 

0 

1379 

3.16 

5.575 

2.25 

1380 

4.3 

4.057 

0 

1381 

4.3 

4.057 

2.25 

1382 

1.93 

4.61 

0 

1383 

1.93 

4.61 

2.25 

1384 

3.31 

3.33 

0 

1385 

3.31 

3.33 

2.25 

1386 

.77 

3.03 

0 

1387 

.77 

3.03 

2.25 

1388 

2.53 

2.356 

0 

1389 

2.53 

2.356 

2.25 

1390 

.17 

1.55 

0 

1391 

.17 

1.55 

2.25 

1392 

2.046 

1.22 

0 

1393 

2.046 

1.22 

2.25 

1394 

0 

0 

0 

1395 

0 

0 

2.25 

1396 

1.87 

0 

0 

1397 

1.87 

0 

2.25 

1398 

8.27 

5.22 

0 

1399 

8.27 

5.22 

2.25 

1400 

8.065 

4.945 

0 

92 


mwm 

1401  8.065  4.945  2.25 

*THE  MESH  FOR  THE  REST  OF  THE  HEADER  IS  GENERATED 


GVOLUME 

1238 

1230 

1231 

1239 

1366 

1367 

1399 

1398 

ELI-2 

EL2-4 

EL3-1 

GVOLUME 

1239 

1231 

1232 

1240 

1398 

1399 

1401 

1400 

ELI -2 

EL2-1 

EL3-1 

GVOLUME 

1240 

1232 

1233 

1241 

1400 

1401 

1369 

1368 

ELI -2 

EL2-1 

EL3-1 

GVOLUME 

1367 

1366 

1370 

1371 

1369 

1368 

1372 

1373 

ELI -2 

EL2-1 

EL3-6 

GVOLUME 

1371 

1370 

1374 

1375 

1373 

1372 

1376 

1377 

ELI-2 

EL2-1 

EL3-6 

GVOLUME 

1375 

1374 

1378 

1379 

1377 

1376 

1380 

1381 

ELI -2 

EL2-1 

EL3-6 

GVOLUME 

1379 

1378 

1382 

1383 

1381 

1380 

1384 

1385 

ELI  *2 

EL2-1 

EL3-6 

GVOLUME 

1383 

1382 

1386 

1387 

1385 

1384 

1388 

1389 

ELI-2 

EL2-1 

EL3-6 

GVOLUME 

1387 

1386 

1390 

1391 

1389 

1388 

1392 

1393 

ELI -2 

EL2-1 

EL3-6 

GVOLUME 

1391 

1390 

1394 

1395 

1393 

1392 

1396 

1397 

ELI-2 

EL2-1 

EL3-6 

’•'THE  NODES  AT  THE  CENTER  OF  THE  TUBES  ARE  FIXED 


BOUNDARIES  111111  TYPE-NODES/ 949  STEP  1  TO  956 
*THE  NODE  AT  THE  ORIGIN  IS  FIXED 
BOUNDARIES  111111  TYPE-NODES / 1892 

*THE  NODES  ABOVE  THE  ORIGIN  ARE  ONLY  ALLOWED  TO  ROLL  IN  THE  X-DIRECTION 
BOUNDARIES  011111  TYPE-NODES 

932  1015  1018  1023  1026  1031  1034  1039  931  1043  93C  1046  929 
1895  1902  1906  1913  1917  1924  1928  1935  1939  1946  1950  1956 
*THE  NODES  IN  THE  Y-0  PLANE  ARE  CONSTRAINED  TO  ROLL  ONLY  IN  THAT  PLANE 
BOUNDARIES  010111  TYPE-NODES 

936  1013  1049  928  961  1051  966  1017  1053  1054  973  1059  978  1025  1061  1062 
985  1067  990  1033  1069  1070  997  1075  935  1041  1077  927  1004  1079  934  1044 
1080  926  1009  1082  933  1047  1083  925  1394  1395  1891  1893  1894  1896  1900 
1901  1903  1904  1905  1907  1911  1912  1914  1915  1916  1918  1922  1923  1925  1926 
1927  1929  1933  1934  1936  1937  1938  1940  1944  1945  1947  1948  1949  1951  1396 
1397  1955  1957 


*THE  NODES  IN  THE  2=0  PLANE  ARE  CONSTRAINED  TO  ROLL  ONLY  IN  THAT  PLANE 
BOUNDARIES  001111  TYPE-NODES 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

17 

28 

39 

50 

61 

72 

83 

94 

105 

116 

127 

133 

136 

139 

142 

145 

148 

151 

157 

163 

166 

169 

175 

181 

184 

187 

190 

193 

199 

205 

208 

211 

217 

223 

226 

229 

232 

235 

241 

247 

250 

253 

259 

265 

268 

271 

274 

277 

280 

283 

286 

289 

292 

295 

301 

307 

310 

313 

319 

325 

328 

331 

334 

337 

343 

349 

352 

355 

361 

367 

370 

373 

376 

379 

385 

391 

394 

397 

403 

409 

412 
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CONSTRAINTS 

*THE  FOLLOWING  CONSTRAINTS  ARE  NECESSARY  BECAUSE  COORDINATES  WERE 

*  DUPLICATED  IN  THE  GENERATION  OF  THE  MESH  FOR  THE  HEADER  AT  THE 

^BOUNDARY  OF  THE  TUBES 

968  1  1019  1 

968  2  1019  2 

968  3  1019  3 

1020  1  1055  1 

1020  2  1055  2 

1056  1  1091  1 

1056  2  1091  2 

1056  3  1091  3 

1127  1  1092  1 

1127  2  1092  2 

1127  3  1092  3 

1159  1  1128  1 

1159  2  1128  2 


967  1  1160  1 
967  2  1160  2 


967  3  1160  3 
980  1  1027  1 


980  2  1027  2 
980  3  1027  3 


1028  1  1063  1 
1028  2  1063  2 


1064  1  1099  1 
1064  2  1099  2 


1064  3  1099  3 
1135  1  1100  1 


1135  2  1100  2 
1135  3  1100  3 


1163  1  1136  1 
1163  2  1136  2 


979  1  1164  1 
979  2  1164  2 
979  3  1164  3 


992  1  1035  1 
992  2  1035  2 
992  3  1035  3 
1036  1  1071  1 
1036  2  1071  2 
1072  1  1107  1 
1072  2  1107  2 
1072  3  1107  3 
1143  1  1108  1 
1143  2  1108  2 
1143  3  1108  3 
1167  1  1144  1 
1167  2  1144  2 
991  1  1168  1 
991  2  1168  2 
991  3  1168  3 
1158  1  1177  1 
1158  2  1177  2 
1158  3  1177  3 
1179  1  1210  1 
1179  2  1210  2 
1125  1  1209  1 
1125  2  1209  2 
1125  3  1209  3 
1211  1  1248  1 
1211  2  1248  2 
1211  3  1248  3 
1284  1  1249  1 
1284  2  1249  2 
1284  3  1249  3 
1320  1  1285  1 
1320  2  1285  2 
1352  1  1321  1 
1352  2  1321  2 
1352  3  1321  3 
1178  1  1353  1 
1178  2  1353  2 
1178  3  1353  3 
1162  1  1185  1 
1162  2  1185  2 
1162  3  1185  3 
1187  1  1216  1 
1187  2  1216  2 
1133  1  1215  1 
1133  2  1215  2 
1133  3  1215  3 
1217  1  1256  1 
1217  2  1256  2 
1217  3  1256  3 
1292  1  1257  1 
1292  2  1257  2 
1292  3  1257  3 
1328  1  1293  1 
1328  2  1293  2 
1356  1  1329  1 
1356  2  1329  2 
1356  3  1329  3 
1186  1  1357  1 
1186  2  1357  2 
1186  3  1357  3 
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1166  1  1193  1 
1166  2  1193  2 
1166  3  1193  3 
1195  1  1222  1 
1195  2  1222  2 
1141  1  1221  1 
1141  2  1221  2 
1141  3  1221  3 
1223  1  1264  1 
1223  2  1264  2 
1223  3  1264  3 
1300  1  1265  1 
1300  2  1265  2 
1300  3  1265  3 
1336  1  1301  1 
1336  2  1301  2 
1360  1  1337  1 
1360  2  1337  2 
1360  3  1337  3 
1194  1  1361  1 
1194  2  ’361  2 
1194  3  1361  3 

*THE  FOLLOWING  CONSTRAINTS 


585  2 

552  2 

585  3 

552  3 

586  2 

553  2 

586  3 

553  3 

587  2 

554  2 

587  3 

554  3 

588  2 

555  2 

588  3 
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589  2 

556  2 

590  2 

557  2 

590  3 

557  3 

591  2 

558  2 
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558  3 

592  2 

559  2 

592  3 

559  3 

593  2 

560  2 

593  3 

560  3 

594  2 

561  2 

594  3 

561  3 

473  2 

470  2 

584  2 

551  2 

584  3 

551  3 

1323  2 

759  2 

1323  3 

759  3 

1354  2 

760  2 

1354  3 

760  3 

1180  2 

761  2 

1180  3 

761  3 

1181  2 

762  2 

1181  3 

762  3 

1182  2 

763  2 

1212  2 

764  2 

1212  3 

764  3 

1213  2 

765  2 

1213  3 

765  3 

1250  2 

766  2 

1250  3 

766  3 

DEFINE  THE  CONTACT  SURFACES 


969  2 

297 

2 

969  3 

297 

3 

970  2 

298 

2 

970  3 

298 

3 

971  2 

299 

2 

971  3 

299 

3 

1021  2 

300 

2 

1021  3 

300 

3 

1022  2 

301 

2 

1057  2 

302 

2 

1057  3 

302 

3 

1058  2 

303 

2 

1058  3 

303 

3 

1093  2 

304 

2 

1093  3 

304 

3 

1094  2 

305 

2 

1094  3 

305 

3 

1129  2 

306 

2 

1129  3 

306 

3 

1130  2 

295 

2 

1161  2 

296 

2 

1161  3 

296 

3 

981  2 

339 

2 

981  3 

339 
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340 

2 
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340 

3 

983  2 

341 

2 

983  3 

341 
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1029  2 

342 

2 

1029  3 

342 

3 

1030  2 

343 

2 

1065  2 

344 

2 

1065  3 

344 

3 

1066  2 

345 

2 

1066  3 
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3 

1101  2 

346 

2 

1101  3 

346 

3 

1102  2 

347 

2 

1102  3 

347 

3 

1137  2 

348 

2 

1137  3 
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3 

1138  2 
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2 

1165  2 

338 

2 

1165  3 
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993  2 
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2 
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1 
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994  3 

995  2 
995  3 

1037  2 

1037  3 

1038  2 
1073  2 

1073  3 

1074  2 
1074  3 
1109  2 

1109  3 

1110  2 
1110  3 
1145  2 
1145  3 
1.146  2 
1169  2 
1169  3 

283  2 

284  2 

284  3 

285  2 

285  3 

286  2 
287  2 

287  3 

288  2 
288  3 

325  2 

326  2 

326  3 

327  2 

327  3 

328  2 

329  2 

329  3 

330  2 
330  3 

367  2 

368  2 

368  3 

369  2 

369  3 

370  2  ] 

371  2  : 

371  3  ] 

372  2  ] 
372  3  J 

409  2  1 

410  2 

410  3 

411  2 

411  3 

412  2  1 

413  2  1 

413  3  1 

414  2  1 
414  3  1 
745  2  1 


3  382  3 

2  383  2 

3  383  3 

2  384  2 

3  384  3 

2  385  2 

2  386  2 

3  386  3 

2  387  2 

3  387  3 

2  388  2 

3  388  3 

2  389  2 

3  389  3 

2  390  2 

3  390  3 

3  379  2 

>  380  2 

!  380  3 

!  1124  2 
!  962  2 

I  962  3 
!  963  2 

I  963  3 

1  1016  2 
1052  2 
1052  3 
1088  2 
1088  3 
1132  2 
974  2 

974  3 

975  2 
975  3 

1024  2 
1060  2 
1060  3 
1096  2 
1096  3 
1140  2 
986  2 

986  3 

987  2 
987  3 

1032  2 
1068  2 
1068  3 
1104  2 
1104  3 
1148  2 
998  2 

998  3 

999  2 
999  3 

1040  2 
1076  2 
1076  3 
1112  2 
1112  3 
1281  2 


747 

2 

1175 

2 

747 

3 

1175 

3 

748 

2 

1176 

2 

749 

2 

1208 

2 

749 

3 

1208 

3 

750 

2 

1245 

2 

750 

3 

1245 

3 

787 

2 

1289 

2 

788 

2 

1325 

2 

788 

3 

1325 

3 

789 

2 

1183 

2 

789 

3 

1183 

3 

790 

2 

1184 

2 

791 

2 

1214 

2 

791 

3 

1214 

3 

792 

2 

1253 

2 

792 

3 

1253 

3 

829 

2 

1297 

2 

830 

2 

1333 

2 

830 

3 

1333 

3 

831 

2 

1191 

2 

831 

3 

1191 

3 

832 

2 

1192 

2 

833 

2 

1220 

2 

833 

3 

1220 

3 

834 

2 

1261 

2 

834 

3 

1261 

3 

871 

2 

1305 

2 

872 

2 

1341 

2 

872 

3 

1341 

3 

873 

2 

1199 

2 

873 

3 

1199 

3 

874 

2 

1200 

2 

875 

2 

1226 

2 

875 

3 

1226 

3 

876 

2 

1269 

2 

876 

3 

1269 

3 

* PORTIONS  OF  THE  MODEL  ARE  DEFINED  FOR  GRAPHICS 

EZONE  TUBEl 

1  2  STEP  1  TO  60 

EZONE  T0BE2 

61  62  STEP  1  TO  120 

EZONE  TUBE3 

121  122  STEP  1  TO  156 

EZONE  ITJBE4 

157  158  STEP  1  TO  192 

ZZONE  COMBO  1  TtJBEl  TUBES 

ZZONE  COMBO 2  TUBE2  TUBE4 

ZZONE  C0MBO3  TUBEl  TUBE2  TUBE3  TUBE4 

EZONE  LOAD 

126  STEP  6  TO  192  203  STEP  12  TO  287  204  STEP  12  TO  288 
ZZONE  LOADING  TUBEl  TUBE2  LOAD 
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